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ABSTRACT 


Simulation  models  were  developed  for  the 
bombardment  of  the  (100),  (110),  and  (111) 
orientations  of  a single-crystal  silicon  lattice  by  an 
argon  ion.  The  interatomic  potentials  included 
repulsive  terms  only,  and  were  approximated  in  closed 
analytic  form  by  least-squares  fitting  of  the 
calculated  interatomic  potentials  obtained  from  the 
Hartree-Foch-Slater  self-consistent  field  equations. 


Sputtering  was  found  to  be  predominantly  a surface 
effect  with  most  sputtering  events  occurring  in  the 
first  three  crystal  layers.  The  (110)  orientation 
demonstrated  a sputtering  mechanism  not  previously 
observed  in  fee  copper  simulations,  and  the  (111) 
orientation  showed  that  reflection  of  the  knock-on 
atom  from  lower  crystal  layers  is  more  common  in 
silicon  than  in  the  fee  copper  lattice.  Adjustment  of 
the  surface  binding  energy  parameter  to  a value  of  4-7 
eV  produced  results  which  agree  with  available 
experimental  data. 
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INTRODUCTION 


The  phenomenon  of  physical  sputtering  has  beer, 
extensively  studied  since  its  discovery  and  description  in 
the  Philosophical  Transactions  by  Grove  in  1852.  Sputtering 
may  be  broadly  viewed  as  the  ejection  of  atomic  material 
from  a target  as  a result  of  bombardment  by  an  energetic 
particle.  The  bombarding  particles  may  be  neutral  atoms, 
ions,  or  the  heavier  elementary  particles;  but  most 
scientific  investigation  nas  been  restricted  to  ionic 
bombardment  because  the  ions  can  be  easily  accelerated  to 
any  desired  velocity.  Sputtering  has  proven  to  be  a 
nuisance  in  the  design  of  high  power  amplifier  and 
oscillator  tubes  because  of  grid  erosion,  in  the  design  of 
thermonuclear  reactors  because  of  plasma  contamination 
caused  by  sputtering  from  the  container  walls,  and  in 
solar-wind  ana  atmospheric  erosion  of  orbiting  satellites. 
Useful  applications  of  sputtering  have  been  found  in  the 
fields  of  ion-getter  pump  design,  thin-film  deposition, 
etching  of  metallurgical  specimens,  and  the  cleaning  and 
etching  of  semiconductor  surfaces.  The  complete  physical 
description  of  the  phenomenon  of  sputtering  is  complicated 
by  limitations  of  experimental  tech  ique,  the  mathematical 
difficulty  of  theoretically  describing  the  behavior  of  a 
many-body  system,  and  the  incomplete  understanding  of  the 
interacting  forces  and  potentials  in  crystal  structures. 

Early  experiments  typically  determined  the  sputtering 
ratio  of  polycrystalline  materials.  The  sputtering  ratio  is 
defined  as  the  number  of  target  atoms  sputtered  per  incident 
ion.  An  exhaustive  study  of  low-energy  sputtering  was 
carried  out  over  the  course  of  seven  years  under  the 
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The  earliest  attempts  to  describe  theoretically  the 
mechanisms  of  sputtering  generally  ignored  the  ordering  of 
atoms  in  the  crystal  lattice.  In  1921,  Thompson  [3] 
proposed  that  atomic  ejection  was  caused  by  the  release  of 
radiation  as  the  ion  struck  the  target.  3ush  and  Smith  [4] 
attempted  to  describe  the  ejection  of  atoms  as  the  result  of 
the  expansion  of  gas  adsorbed  by  the  target  material.  Von 
Hippel  [5]  presented  the  'Hot-Spot  Theory'  in  1926,  which 
stated  that  the  energy  dissipation  of  the  slowing  ion 
evaporated  target  material  from  the  surface.  A theory  based 
on  momentum  transfer  between  the  bombarding  ion  and  the 
lattice  atoms  was  proposed  by  Lamar  and  Compton  [5]  in  1934. 
Their  theory  stated  than  ions  may  penetrate  the  material  and 
recoil  from  lower  layers,  causing  ejection  of  surface  atoms. 
In  1956  Harrison  [7]  proposed  a theory  for  sputtering  from 
amorphous  materials  based  on  transport  theory.  Silsbee  [9] 
proposed  the  theory  that  sputtering  in  preferred  directions 
from  single  crystals  is  the  result  of  the  focusing  of 
momentum  along  a line  of  atoms  in  the  crystal.  This  theory 
was  questioned  by  Lehmann  and  Sigmund  [91  because  of  the 
extremely  short  pathlengt'n  of  the  Silsbee  focusons  even  at 
high  energies.  In  1969,  Sigmund  [10]  reported  his  theory 
for  the  sputtering  of  amorphous  and  polycrystalline 
materials.  His  theory  was  also  based  on  the  Boltzmann 
transport  equation. 
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Analytic  sputtering  theories  cannot,  in  general, 
adequately  define  and  treat  the  parameters  and  variables 
required  in  the  theory  of  single-crystal  sputtering.  The 
development  of  high-speed  computers  offered  the  chance  to 
better  understand  the  sputtering  process  through  the  use  of 
simulation  models.  In  1960,  Gibson,  Goland,  Milgram,  and 
Vineyard  [11]  built  a computer  model  to  represent  metallic 
copper,  and  studied  radiation  damage  events  at  energies  up 
to  400  eV.  In  this  model,  one  atom  in  the  stationary 
lattice  was  given  an  arbitrary  kinetic  energy  and  direction 
of  motion,  and  was  allowed  to  interact  with  the  atoms  in  the 
lattice,  resulting  in  a cascade  of  independent  binary 
collisions.  In  1967,  Harrison,  Levy,  Johnson,  and  Effron 
[12]  used  a computer  to  simulate  the  bombardment  of  a single 
fee  copper  crystal  by  argon.  Through  computer  simulation, 
the  sputtering  mechanisms  in  the  fee  copper  crystal  were 
identified.  Harrison  et  al  reported  that  for  ion  energies 
less  than  10  keV,  the  sputtering  process  in  fee  crystals 
occurred  predominantly  within  three  atomic  layers  of  the 
surface.  The  initial  computer  simulation  included  only  the 
repulsive  forces  between  the  atoms.  The  algorithm  utilized 
to  solve  the  equations  of  motion  in  this  simulation  and 
later  simulations  was  described  in  1969  by  Harrison,  Gay, 
and  Effron  [13].  Continuing  efforts  to  develop  a more 
precise  computer  model  resulted  in  the  development  of  a 
Cu-Cu  potential  function  which  included  an  attractive 
portion,  inclusion  of  surface  layer  relaxation  in  the 
crystal  model,  and  determination  of  surface  atom  binding 
energies  of  the  (100),  (110),  and  (111)  Cu  crystal 
orientations.  In  1972,  Harrison,  Uoore,  and  Holcombe  [14] 
showed  improved  agreement  between  simulation  results  and 
experimental  data,  but  the  original  interpretation  of  the 
low-energy  sputtering  mechanisms  in  fee  crystals  remained 
unchanged. 
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This  thesis  is  an  extension  of  the  previous  NPS  computer 
simulation  studies  to  silicon  bombarded  by  argon.  Silicon, 
a group  IV  element,  forms  a crystal  in  the  diamond  lattice 
structure.  An  acceptable  computer  model  of  single-crystal 
silicon  bombardment  must  generate  the  lattice  in  the  proper 
orientation,  and  it  must  simulate  with  reasonable 

approximations  of  the  atomic  interaction  potentials  and 
forces  for  both  the  Ar-Si  and  Si-Si  collisions.  A 

statistically  representative  impact  area  of  the  Si  lattice 
must  be  determined  for  each  crystal  orientation,  a suitable 
force  and  potential  truncation  distance  (which 
reasonable  computer  run  time  while  ensuring 
principle  of  energy  conservation  is  not  violated) 
determined,  and  the  maximum  allowable  distance  m: 
most  energetic  atom  in  any  given  timestep  must  be 
so  that  forces  and  velocities  are  updated  often 
ensure  energy  conservation  without  ma 
prohibitively  long.  For  this  prel 

model,  only  the  repulsive  portion 
potentials  is  included,  and  surface 
considered.  Furthermore,  the  lack  of 
yields  of  single-crystal  silicon  leaves  the 
surface  binding  energy  unresolved, 
required  to  free  an  atom  from  the 
binding  energy.  For  sputtering  to  oc 
of  a lattice  atom  must  exceed  the 
and  the  atom  must  be  moving  away  from 
Adjustment  of  this  parameter  to 
experimental  data  results  in  a m 
However,  this  parameter  has  no  effect 
mechanisms  involved  in  the  model, 
whenever  suitable  data  become  available 
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The  diamond  lattice  model  developed  possesses  the  same 
features  as  the  original  fee  copper  model  developed  by 
Harrison  et  al  [12].  Further  refinements  of  this  model 
should  include  surface  relaxation,  interatomic  cohesive 
energy,  and  perhaps  non-central  force  and  potential 
approximations.  From  earlier  experience  with  the  copper 
crystal  however,  it  is  expected  that  this  preliminary 
identification  of  sputtering  and  momentum  transfer 
mechanisms  will  not  be  affected  greatly  by  the 
overs  imp lif ications. 


II. 


MODEL  DEVELOPMENT 


A.  GENERATION  OF  THE  SILICON  DIAMOND  LATTICE 

The  directional  covalent  bonds  between  a silicon  atom 
and  its  four  nearest  neighbors  result  in  the  formation  of  a 
symmetrical  tetrahedron  (Fig.1).  The  bonding  of  these 
tetrahedra  with  successive  nearest  neighbors  then  results  in 
the  formation  of  the  diamond  lattice  (Fig  2) . The  lattice 
constant  a in  silicon  is  5.43  angstroms,  and  the  lattice 
unit,  which  is  the  width  of  the  basic  tetrahedron,  is  2.72 
angstroms.  In  the  development  of  the  computer  model,  the 
lattice  unit  ( LU ) proves  to  be  a convenient  unit  of  length, 
and  is  used  extensively.  This  lattice  is  fee  with  a basis 
of  two  atoms  associated  with  each  lattice  point.  The  second 
atom  is  displaced  0.86602  LU  along  the  (111)  diagonal  from 
the  first.  Figure  3 shows  the  arrangement  of  the  silicon 
lattice  in  the  (100)  orientation.  Square  abed  in  Figur°  3 
encloses  a statistically  representative  impact  area  for 
normally  incident  ions,  which  by  folding  along  two  axes  of 
symmetry  (e.g.  ab  and  ad)  will  map  into  all  possible  impact 
points  on  the  (100)  surface.  The  (100)  channel  is  shown  at 
the  center  of  square  abed.  The  distance  from  the  center  of 
this  channel  to  the  center  of  the  nearest  atom  in  the  same 
vertical  plane  is  0.5  LU.  Figure  4 shows  the  specific 
impact  points  chosen  for  the  (100)  surface. 

The  (110)  orientation  of  the  diamond  lattice  reveals  a 
large  channel  (Fig.  5)  through  which  the  incoming  ion  may 
pass  without  making  a hard  collision  with  any  lattice  atom. 
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The  distance  from  the  center  of  this  channel  to  the  center 
of  the  closest  atom  in  the  same  plane  is  0.75  LU.  This 
distance  is  the  maximum  impact  parameter  seen  in  any  of  the 
three  orientations,  and  the  Si-Ar  potential  and  force 
functions  must  be  modeled  at  least  to  this  distance.  Figure 
6 shows  the  representative  impact  area  used  for  normally 
incident  ions  on  the  (110)  plane.  The  impact  area  abed 
contains  30  Impact  points,  and  all  possible  impact  points 
may  be  represented  by  a folding  about  the  axis  of  symmetry 
ab. 


A further  rotation  of  the  lattice  yields  the  (ill) 
orientation  (Fig.  7).  The  large  channel  evident  in  the 
(110)  orientation  is  no  longer  present,  but  an  ion  normally 
incident  at  point  b may  travel  relatively  deep  into  the 
lattice  (the  fifth  plane)  before  experiencing  more  than  a 
grazing  collision.  The  small  (111)  channels  surrounding  the 
fifth  plane  atom  have  an  impact  parameter  of  0.471  LU.  This 
may  be  viewed  as  a "pseudo-channeling"  effect,  and  suggests 
that  perhaps  the  (110)  and  (111)  surfaces  may  exhibit 
somewhat  similar  behavior.  A more  complicated  folding  of 
any  triangle  in  the  hexagon  will  yield  the  total  area 
abedefg.  Triangle  bed  was  chosen  for  the  computer  model. 
Figure  3 shows  the  position  of  the  36  impact  points  used  in 
the  (111)  orientation. 

The  subroutines  utilized  to  generate  the  three  lattice 
orientations  are  a modification  of  the  lattice  generators 
used  by  Finno  in  silicon  channeling  studies  [15],  and  are 
discussed  in  Appendix  A. 
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B.  INTERACTION  POTENTIALS  AND  FORCE  FUNCTIONS 


In  1967,  Wedepohl  *16]  reported  the  calculation  of  the 
interaction  energies  between  like  atoms  and  ions  in  the 
energy  range  which  is  of  interest  in  radiation  damaae 
calculations  (i.e.,  50  e?-105  eV) . Wedepohl  compared  the 
interaction  energies  obtained  from  three  theoretical  models: 

(i)  the  guantu m- mechanical  (QM)  radial  electron  density 
obtained  using  Hartree  self-consistent  field  calculations 
(see  ref . 17), 

(ii)  Thomas-Fer mi  (TF)  electron  distributions  (see  ref. 

18)  . 

(iii)  Thomas-Fermi-Dirac  (TFD)  electron  distributions 
(see  ref . 19). 

An  unpublished  modification  of  the  Wedepohl  QM  method 
developed  by  Harrison  was  utilized  to  obtain  the  Si-Si  and 
Ar-Si  interaction  potential  functions.  Tabulated  values  of 
the  radial  electron  densities  of  various  ionic  states 
(including  the  neutral  states)  of  Ar  and  Si  were  obtained 
from  the  Her man-Skillman  computer  program  [20],  which  is 
based  on  the  Slater  approximation  of  the  Hartree-Foch 
self-consistent  field  calculations.  These  densities  were 
then  used  in  a program  similar  to  Wedepohl's  [16]  to  obtain 
tabulated  values  of  the  interaction  potentials.  It  should 
be  kept  in  mind  that  the  adiabatic  approximation  was 
utilized,  i.e.  the  electronic  charge  distributions  were  not 
allowed  to  deform  due  to  electronic  repulsion  as  the  atoms 
were  brought  in  close  proximity  to  each  other.  This  results 
in  a screening  function  which  is  in  all  likelihood  too  high, 
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yielding  a higher  (harder)  repulsive  potential  between  atoms 
than  would  be  obtained  if  the  electron  clouds  were  allowed 
to  deform.  In  addition,  since  a radial  density  function  was 
utilized  and  angular  dependence  was  ignored,  the  resulting 
potential  does  not  include  the  directional  properties  of  the 
Si  electronic  wave  functions.  Also,  since  only  the 
repulsive  portion  of  the  potential  is  calculated,  the  Si-Si 
potential  must  be  truncated  at  the  lattice  nearest-neighbor 
distance  in  order  to  prevent  spontaneous  crystal  expansion 
and  disintegration.  As  was  stated  earlier,  the  truncation 
point  of  the  Ar-Si  potential  must  be  at  a distance  greater 
than  0.75  LU.  The  truncation  distance  was  chosen  to  be  the 
same  as  the  nearest- neighbor  distance  (0.86602  LU)  . 

Figures  9 and  10  show  the  Si  and  Ar  radial  electron 
densities  obtained  from  the  Her.' nan-S ici lima n program.  Figure 
11  shows  the  potentials  obtained  from  the  interaction  of  Ar 
and  Si  + 4,  and  from  the  interaction  of  Ar+1-Si.  The 
Ar-si+4  interaction  was  chosen  for  the  model  since  it  is 
lower  at  the  nearest-neighbor  distance.  It  is  not 
unreasonable  to  assume  that  under  experimental  conditions 
many  of  the  accelerated  Ar  ions  will  be  neutralized  from  the 
conduction  electrons  of  the  metal  before  collision  with  the 
first  lattice  atom.  The  silicon  atom  in  the  lattice  will 
appear  as  an  ion  to  the  approaching  argon  atom  because  of 
the  covalent  sharing  of  its  outer  shell  electrons.  Since 
the  potential  function  in  all  likelihood  is  too  nard,  the 
choice  cf  a lower  potential  is  reasonable. 

Figure  12  shows  the  selected  Si-Si  potential.  The 
choice  of  Si-Si+l  was  made  in  order  to  obtain  a potential 
function  which  closely  approached  zero  at  the  equilibrium 
separation  (nearest-neighbor  distance)  in1  the  lattice. 

The  potential  values  obtained  were  then  approximated  in 
a closed  analytical  form  in  order  that  they  could  be  easily 
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utilized  in  the  sputtering  program.  Wilson,  Haggmark,  and 
Biersack  [21  ] recently  reported  the  approximation  of  14 
diatomic  interactions  by  the  least-sguares  fitting  of  the 
f ree-electron  potential  solutions  with  a Moliere-like 
potential  form: 

3 

V(r)  = (Z,Z2e2A)  I C;  exp(-bif/a) 

1 * i =i 

where  (a)  is  a screening  length  and, 

3 

EC;  = I 

i = l 

This  method  averaged  approximately  10%  accuracy,  and  is  seen 
to  resemble  the  method  used  in  the  development  of  this 
sputtering  model. 


1.  Ar-Si+4  Interaction  Potential 


Examination  of  the  Ar-Si*4  potential  plot  showed 
that  for  a separation  distance  greater  than  0.25  LU,  the 
logrithmic  plot  was  essentially  a straight  line  suggesting  a 
Born-Mayer  potential  function  of  the  form: 

V (r)  = ex  p[  PE  XA  ♦ PEX3*r]. 

3M 

A least  squares  fit  of  the  region  between  0.25  LU  and  0.75 
LU  (the  maximum  impact  parameter)  yielded  a value  for  these 
exponential  coefficients: 

V (r)  = exp[  8. 40056-10. 85482r], 

BM 

where  V (r)  is  potential  in  electron  volts  and  r is 
separation  distance  in  lattice  units.  A "peeling"  process 
was  then  utilized,  where  the  value  of  the  Born-Mayer  portion 
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of  the  potential  was  subtracted  from  each  data  point.  This 
process  revealed  that  the  region  of  separation  less  than 
0.10916  LU  could  be  represented  by  a Bohr  potential  of  the 
form: 


V (r) -V  (r)  = (1/r)  exp[  A + B«r]. 

B M 

A least  squares  fit  of  the  region  less  than  0.10916  LU 
yielded  the  following  exponential  coefficients: 

V (r )-V  (r)  = (1/r) exp[  7.  164  10-2  1 . 10187r  ]. 

BM 

In  order  to  insure  continuity  of  the  analytic  form  of  the 
potential  function  and  its  first  derivative  with  respect  to 
r,  a "transition"  function  was  utilized  at  the  cross-over 
point  (0.10916  LU) . This  transition  function  proved  to  be 
of  the  Born-Mayer  form  also: 

V = exp[  10.36930-30.  13574r  ]. 

D 

The  modeling  of  the  Ar- Si+4  potential  function  was 
complete,  yielding  a function  of  the  form: 


m/r)  «xp  (7164IO-2I  IOI87r)  ♦ VBM(r);  r<OIOI96  LU 
v(r)  = 1 exp  (10  36930-3  0 13574  r)  ♦ VBM(r)i  r >0*10196  LU 

Figure  13  is  the  plot  of  the  analytic  model  of  the  potential 
function  and  corresponding  potential  values  obtained  from 
the  Herman-Skiilraan  and  modified  Wadepohl  programs.  The 
maximum  relative  error  in  the  analytic  potential 

representation  between  0.0  and  0.75  LU  occurs  at  0.24950  LU, 
where  the  analytic  model  is  1.67%  higher  than  the 
corresponding  computed  value.  This  is  caused  by  the 

"transition  potential"  which  has  not  damped  out  sufficiently 
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at  this  point,  but  which  rapidly  disappears  into  the 
Born-flayer  potential.  The  potential  model  obtained  was 
incorporated  into  the  ENERGY  subroutine  of  the  basic 
computer  program  (Appendix  B)  as  a two-branch  potential 
function  with  cross-over  point  at  0.10916  LU  and  truncation 
at  0.86602  LU  (nearest-neighbor  distance). 


2.  Si-S i+ 1 Potential  Function 


Examination  of  the  Si-Si+1  potential  values  showed 
that  for  separation  distances  above  0.5458  LU  the  potential 
function  was  essentially  of  phe  Born-flayer  form: 

V (r)  = expfSXA  + EXE*r  ]. 

BM 

A two  v aria  tie  linear  regression  yielded  the  following 
coefficients: 

V (r)  = exp[  6.  68394-6 .94475r], 


where  V (r)  is  in  eV,  and  r is  in  LU.  "Peeling"  of  this 
BM 

potential  from  the  tabulated  values  yielded: 


V (r)  -V  (r)  = (1/r)exp[  6. 95261-  15. 56752r] 
BM 


Continuity  was  obtained  without  the  use  of  a "transition 
potential"  by  placing  the  crossover  point  at  0.7  LU. 
Figure  14  shows  the  analytic  form  of  the  Si-Si+l  potential 
and  corresponding  Herman-skillraan  data  points.  The  relative 
error  in  the  analytic  form  of  the  potential  function  is  less 
than  six  percent  for  all  energies  less  than  10.0  keV. 


l 
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This  potential  representation  was  also  incorporated  into  the 
ENERGY  subroutine  (Appendix  3),  with  truncation  at  the 
nearest-neighbor  distance. 

3 • Ar-Si+4  and  S i -S i+i  Force  Functions 

The  force  functions  were  simply  determined  by 
obtaining  the  negative  gradients  of  the  potential  functions: 

F(r)  = - a V(r)/ dr 

The  problem  of  an  infinite  gradient  at  the  potential 
truncation  point  was  handled  as  demonstrated  by  Harrison  et 
al  [13].  By  utilizing  the  average  force: 

< F;>  = [ V(x)-  V(R)]/(R-x)-,  R-D<x<  R 

where  E is  the  nearest-neighbor  distance  and  D is  the 
maximum  allowed  displacement  of  the  most  energetic  atom  in 
one  timestep,  the  force  is  gently  "turned  on"  at  tne 
nearest-neighbor  distance.  The  value  for  D which  allows 
both  reasonable  run  times  and  small  energy  decrement  was 
found  to  be  0.075  LU.  The  two-branch  force  approximations 
were  incorporated  into  the  program  in  subroutine  STEP 
(Appendix  C) . 


L 
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The  three  lattice  orientation  models  were  utilized  in 
computer  program  runs  at  600-eV  and  1-keV  ion  bombardment 
energies.  The  (110)  orientation  model  was  also  utilized  in 
2-keV  ion  bombardment  programs.  The  program  runs  at  all 
energies  revealed  no  large  energy  decrements  from  timestep 
to  timestep,  leading  to  the  conclusion  that  the  choice  of 
0.075  LU  for  D was  reasonable.  A computer  run  with  D set  at 
0.1  LU  showed  a decreased  sputtering  ratio  and  a degradation 
of  the  ability  to  maintain  energy  conservation.  The 
preliminary  investigative  runs  of  the  silicon  model  were 
made  utilizing  a crystal  size  of  8x4x3  in  the  case  of  the 
(100)  and  (111)  orientations,  and  8x8x8  in  the  case  of  the 
(110)  orientation.  Each  model  contained  256  atoms  (the 
(110)  orientation  generates  only  1 atom  per  lattice  site 
while  the  (100)  and  (111)  orientation  generate  two,  hence 
the  apparant,  but  non-real,  increase  in  the  lattice  size). 
The  mean  elapsed  computer  time  required  to  run  a completely 
compiled  and  linked  load  module  on  the  Naval  Postgraduate 
School's  IBM  360/67  computer  was  approximately  160  minutes. 
The  lattice  size  of  256  atoms  did  not  sufficiently  contain 
all  the  possible  events  even  at  these  low  ion  bombardment 
energies;  however,  the  size  of  the  lattice  was  deemed 
sufficient  tc  test  the  model  while  still  minimizing  computer 
run  times.  The  sputtering  ratios  were  determined  for 
lattice  binding  energies  of  1 eV  to  10  eV  in  1-eV 
increments,  and  also  at  15  eV  for  each  orientation.  The 
mechanisms  of  sputtering  were  then  traced  at  active  impact 
points  in  each  orientation. 


Prior  to  the  development  of  the  diamond  lattice  model, 
work  had  been  done  on  a fee  model  in  which  the  size  of  the 
lattice  expanded  when  atoms  near  the  boundaries  began  to 
move.  The  model  was  developed  but  was  unsuccessful  because 
momentum  transfer  to  the  boundary  was  so  rapid  that  the 
crystal  expanded  to  maximum  size  very  shortly  after  ion 
impact,  and  no  computer  time  was  saved  in  the  simulation. 
The  spread  of  momentum  in  a fee  crystal  was  compared  with 
the  spread  in  the  diamond  lattice  to  determine  the 
differences  in  the  momentum  spreading  mechanism,  and  to 
determine  the  feasibility  of  an  expanding  diamond  lattice 
model.  The  momentum  spreading  mechanism  was  found  to  be 
different,  but  the  rapid  spread  of  momentum  to  the 
boundaries  also  was  observed  in  the  diamond  lattice. 

Although  the  investigation  of  the  ion  bombardment  of 
silicon  is  by  no  means  complete,  the  results  obtained 
indicate  that  the  model  is  physically  reasonable,  and 
further  investigation  coupled  with  correlation  with 
experimental  data  obtained  from  other  sources  could  prove 
fruitful  in  enhancing  the  understanding  of  the  problems 
involved  in  the  etching  and  ion  implantation  of  silicon  and 
germanium  semiconductor  materials. 

A.  (100)  LATTICE  ORIENTATION 

Figure  15  shows  a roughly  isometric  projection  of  the 
256  atoms  in  an  8x4x8  diamond  lattice  with  (100) 
orientation.  This  type  of  projection,  first  used  by 
Harrison  and  Delaplain  [22],  with  the  distance  between 
layers  in  the  Y direction  exaggerated,  and  the  atoms 
represented  by  20°  ellipses  has  proved  to  be  a valuable  tool 
for  analyzing  the  computer  programs  and  for  displaying  data. 
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The  first  atom  in  the  lattice  is  atom  number  three,  with 
atom  number  one  being  the  Ar  ion,  and  atom  number  two 
reserved  for  other  projects.  Future  studies  may  utilize 
atom  number  two  as  an  interstitial  atom,  surface 
irregularity  (adatom) , or  the  second  atom  in  a diatomic 
bombarding  molecule  without  disturbing  the  numbering  of  the 
basic  crystal.  In  the  display  of  sputtered  atoms  and  in  the 
tracing  of  ejection  mechanisms,  lower  planes  which  were  not 
directly  involved  in  the  event  were  eliminated  from  the 
figures,  but  included  in  all  calculations. 

1.  600-eV  Results 


The  complete  4 1 - impact-point , 600-eV  simulation 
utilized  175  minutes  of  computer  time  on  an  8x4x8  size 
crystal.  Kith  the  surface  binding  energy  (ETHR)  set  at  1 
eV,  123  atoms  were  sputtered  yielding  a sputtering  ratio  of 
3.00.  Figure  16  shows  that  the  low  energy  sputtering  of  the 
(100)  orientation  is  basically  a surface  event,  with  46  of 
the  atoms  in  the  first  layer  being  sputtered,  11  atoms  in 
the  second  layer  being  sputtered,  and  only  one  sputtering 
event  (atom  149)  taking  place  in  the  third  layer.  Atom  149 
was  ejected  with  a velocity  normal  to  the  crystal  face 
corresponding  to  a kinetic  energy  of  1.33  eV,  a low  energy 
ejection.  All  atoms  that  were  sputtered  at  more  than  three 
impact  points  were  contained  in  the  first  layer,  with  atom 
32  being  sputtered  seven  times,  atom  39  being  sputtered 
eight  times,  atom  40  being  sputtered  five  times,  and  atom  55 
being  sputtered  five  times.  In  the  lower  layers,  atom  95 
was  sputtered  three  times,  and  atom  103  was  sputtered  twice. 
All  other  sputtered  atoms  in  the  lower  layers  sputtered  only 
once.  Atom  103  is  inside  the  impact  area,  and  atom  95  is 
directly  adjacent  to  it;  thereby  receiving  large  kinetic 
energies  very  soon  after  impact  of  the  ion  the  lattice. 
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A plot  of  the  sputtering  ratio  vs.  the  surface 
binding  energy  was  developed  in  order  to  gain  insight  into 
the  value  of  the  surface  binding  energy.  Figure  17  shows 
the  results  of  this  plot  for  the  600  eV  program.  The  rapid 
increase  in  sputtering  ratio  below  3 eV  indicates  that  the 
surface  binding  energy  is  higher  than  this  value.  The  more 
gently  sloping  decrease  in  sputtering  ratio  above  3 eV 
yields  little  insight  into  the  value  of  the  binding  energy. 
The  sputtering  ratio  of  0.122  at  a binding  energy  of  15  eV 
indicates  that  5 atoms  were  ejected  with  a normal  velocity 
component  greater  than  13,152  m /sec.  The  most  energetic 
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energy  greater  than  15  eV  occurea  in  the  first  and  second 
layers  (four  in  the  first  layer,  one  in  the  second  layer). 

The  most  active  impact  point  in  the  600-e/  program 
occurea  when  atoms  39  and  40  were  impacted  at  point  5025 
(impact  point  displaced  0.50  LU  in  X direction  from  the 
center  of  atom  39,  and  0.25  LU  in  Z direction,  hence  5025). 
Eight  atoms  were  sputtered  at  this  point.  A single- impact 
program  was  run  at  this  point  with  a listing  of  the 
positions,  velocities,  kinetic  and  potential  energies,  and 
forces  exerted  on  each  atom  which  possessed  a total  energy 
greater  than  0.1  eV.  The  listing  was  printed  once  every 
five  timesteps.  The  atom  motions  which  caused  the  eight 
atoms  to  be  sputtered  were  traced  in  order  to  gain  insight 
into  the  sputtering  mechanisms  involved  in  the  (100)  diamond 
lattice  orientation.  The  advantage  of  computer  simulation 
in  which  each  atom  is  uniquely  identified  is  evident  in  the 
following  mechanism  summary. 
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Mechanism  Summary 


Atoms  7,9,15,18,32,57,65,  and  119  were  sputtered 
as  a result  of  an  ion  impact  at  impact  point  5025.  Atom  57 
was  sputtered  at  timestep  100,  atom  65  at  timestep  105, 
atoms  7,9,  and  15  at  timestep  110,  atom  18  at  timestep  115, 
and  atom  32  was  sputtered  at  timestep  130.  Each  sputtered 
atom  was  followed  separately,  and  the  individual  traces  are 
illustrated  in  Figure  18. 

(1)  Atom  32.  The  impacting  Ar  ion  (atom  1) 
traveled  into  the  second  crystal  layer  after  striking  atom 
40.  At  timestep  20,  atom  1 passed  within  0.508  L(J  of  atom 
103,  causing  103  to  move  in  the  positive  X direction.  As 
atom  103  moved  below  atom  32,  it  forced  32  upward.  Atoms  31 
and  41  moved  upward  and  away  from  atom  32  as  it  was  still 
gaining  velocity  from  atom  103.  At  timestep  45,  atom  103 
reached  its  point  of  closest  approach  with  atom  32.  As  a 
result  of  the  32-103  collision,  atom  32  gained  an  upward 
velocity  (having  been  slowed  slightly  by  grazing  collisions 
with  31  and  41)  . At  timestep  60  atoms  32  and  103  were  no 
longer  interacting,  and  atom  32  had  risen  0.118  LU.  Atom  32 
continued  to  rise,  interacting  with  no  other  atoms  until 
timestep  125.  At  that  time  atom  32  interacted  with  atom  25 
which  was  rising  with  a low  velocity  as  a result  of  atom 
movement  below  it.  This  slight  interaction  slowed  atom  25, 
but  raised  atom  32's  normal  velocity  component.  At  timestep 
130,  atoms  25  and  32  were  no  longer  interacting,  and  atom  32 
continued  to  rise  with  no  further  interactions.  The 
mechanism  observed  in  the  sputtering  of  atom  32  was  the 
grazing  collision  with  a lower  atom  in  the  basic  tetrahedron 
composed  of  atoms  31,32,41,97,  and  103.  The  velocity 
imparted  to  atom  103  (roughly  parallel  to  the  X,Z  plane) 
resulted  in  the  rising  of  the  center  and  upper  atoms  in  the 
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tetrahedron . 
a high  enough 
sput.  This 
bombardments . 


Atom  32,  however,  was  the  only  atom  possessing 
normal  velocity  component  to  be  counted  as  a 
may  not  always  be  the  case  in  higher  energy 


(2)  Atom  _15.  Atom  1 in  its  passage  through  the 
second  layer  also  interacted  with  atom  96,  causing  it  to 
move  in  the  Z direction  with  a velocity  component  of  -29,516 
m /sec  while  traveling  downward  at  6715  m/sec.  The  atoms 
below  96  slowed  its  downward  velocity  slightly.  At  timestep 
60  atom  80  was  struck  by  atom  96,  which  caused  atom  80  to 
rise  rapidly.  As  96  continued  to  pass  under  80,  it  imparted 
a high  negative  Z velocity  component  to  80.  Atom  80  then 
rose  under  atom  15.  Atom  15  interacted  with  no  other  atom, 
quickly  gained  a Y velocity  component  of  -7838  m/sec  and  was 
counted  as  sputtered  at  timestep  110.  Here  a mechanism 
similar  to  the  scoop  or  squeeze  mechanism  in  fee  copper 
described  by  Harrison  et  al  [12]  took  place  in  the  second 
layer,  causing  a surface  atom  to  be  sputtered. 

(3)  Atom  5.  After  collision  with  atom  15,  atom 
80  was  deflected  below  atom  9,  reaching  the  point  of  closest 
approach  at  timestep  100.  By  timestep  105,  atoms  80  and  9 
were  no  longer  interacting  and  atom  9 rose  from  the  surface 
with  a normal  velocity  component  of  -9311  m/sec,  interacting 
with  no  other  atom.  At  timestep  110,  atom  9 was  also 
counted  as  a sput.  This  event  was  also  a result  of  the 
scoop  or  squeeze  mechanism  in  the  second  layer.  Atom  80  was 
deflected  upward  by  atom  10,  and  also  rose  from  the  surface 
with  a velocity  slightly  too  small  to  be  considered  a 
sputtered  atom.  This  illustrates  a case  of  a second  layer 
atom  rising  from  the  surface  through  a vacancy  created  by 
impact  from  lower  layers,  and  it  is  reasonable  to  assume 
that  at  other  times  this  occurance  causes  sputtering. 
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(4)  Atom  2*  As  atom  80  rose  under  atom  15,  it 
also  interacted  with  atom  8,  reaching  the  point  of  closest 
approach  at  timestep  80.  Atom  8 traveled  in  the  negative  X 
and  Z directions  while  rising  into  collision  with  atom  7. 
The  motion  of  atom  8,  which  was  approximately  at  an  angle  of 
45°  with  the  normal  to  the  surface,  caused  atom  7 to  be 
accelerated  in  the  same  direction.  At  timestep  100,  atom  7 
no  longer  interacted  with  the  greatly  slowed  atom  8,  but 
rose  from  the  surface  with  an  X velocity  of  -4893  m/sec,  a Y 
velocity  of  -4559  m/s ec,  and  a Z velocity  of  -4338  ra/sec. 
At  timestep  110,  atom  7 was  also  counted  as  a sput.  The 
scoop  or  squeeze  mechanism  in  the  second  layer  caused  the 
center  atom  in  tetrahedron  7,8,15,71,79  to  impart  momentum 
to  its  upper  atom  at  an  angle  similar  to  the  deep  mechanism 
described  by  Harrison  et  al  [12];  however,  the  mechanism  is 
unique  to  the  diamond  lattice  where  the  spread  of  momentum 
through  the  tetrahedral  atoms  occurs  at  all  depths,  and 
seemingly  in  all  directions. 

(5)  Atom  1.8.  After  traveling  under  atom  9, 
atom  80  collided  with  atom  10  at  timestep  80,  and  then  rose 
upward  from  the  surface.  Atom  10  passed  under  atom  17, 
causing  it  to  rise  at  a velocity  too  low  to  be  considered  a 
sputtered  atom,  and  then  passed  below  atom  18  at  a distance 
of  0.803  LU,  imparting  a velocity  of  1891  m/sec  in  the  X 
direction,  -2359  m/sec  in  the  Y direction,  and  3389  m/sec  in 
the  Z direction  to  atom  13.  Atom  18  then  crossed  the 
vacancy  between  it  and  atom  34,  colliding  with  atom  34  at 
timestep  150.  As  atom  18  passed  above  34,  its  normal 
velocity  component  was  increased  to  -3689  m/sec,  while  its  Z 
component  was  reversed  to  375  m/sec.  At  timestep  155,  atom 
18  was  counted  as  a sputtered  atom,  also  having  been 
sputtered  by  the  second-layer  mechanism  affecting  the  motion 
of  atom  80. 
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(6)  Atom  U.2*  Atom  1 struck  atom  40  as  it 
entered  the  crystal,  driving  atom  40  downward  into  collision 
with  atoms  112  and  119.  Atom  119  was  then  driven  toward 
atom  120,  and  at  timestep  50  was  driven  upward  by  the 
combination  cf  interactions  with  atoms  40  and  120.  At 
timestep  70,  atom  119  was  no  longer  interacting  with  atoms 
40  and  120,  but  was  deflected  upward  further  by  atom  129, 
and  was  deflected  in  the  negative  X and  positive  Z 
directions  by  collision  with  atoms  57  and  58.  At  timestep 
90,  atom  119  was  interacting  with  no  other  atoms,  and  was 
traveling  with  an  X velocity  of  -2149  m/sec,  a Y velocity  of 
-4215  m/sec,  and  a Z velocity  of  3460  m/soc.  At  timestep 
105,  atom  119  collided  with  atom  64,  reversing  its  X 
velocity  and  slowing  its  negative  Y velocity.  3y  timestep 
115,  atom  119  had  risen  to  0.157  L0  above  the  surface,  and 
was  interacting  with  no  other  atom.  At  timestep  130,  atom 
119  was  counted  as  a sputtered  atom,  possessing  a normal 
velocity  component  of  -3035  m/sec.  Again  the  squeeze 
mechanism  caused  an  atom  to  rise  from  the  second  layer  and, 
in  this  case,  move  into  the  (100)  channel  formed  by  atoms 
57,58,119,  and  120,  be  reflected  out  of  the  the  channel  and 
rise  due  to  a secondary  collision  with  another  surface  atom. 

(7)  Atom  57.  Atom  119  while  rising  to  the 
surface,  collided  with  atoms  57  and  58  at  timestep  70.  At 
timestep  75,  atom  119  passed  within  0.762  LU  of  atom  57 
causing  it  to  rise  from  the  surface.  At  timestep  90,  atom 
57  was  no  longer  interacting  with  atom  119,  and  possessed  a 
normal  velocity  component  of  -4312  m/sec.  Atom  57  continued 
to  rise,  interacting  with  no  other  atoms,  and  was  counted  as 
a sputtered  atom  at  timestep  100.  Again  the  squeeze 
mechanism  in  the  second  layer  resulted  in  the  sputtering  of 
a surface  atom. 
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(8)  Atom  §5.  The  collision  of  atom  119  with 
atom  58  resulted  in  the  passing  of  atom  59  below  atom  65  at 
a distance  of  0.748  LU  at  timestep  85.  5y  timestep  95,  atom 
65  was  no  longer  interacting  with  atom  58,  and  continued  to 
rise  with  no  further  interactions  until  being  counted  as  a 
sputtered  atom  at  timestep  105.  Atom  65  possessed  a normal 
velocity  component  of  -4185  m/sec  after  interaction  with 
atom  58.  Once  again  the  squeeze  mechanism  in  the  second 
layer,  combined  with  momentum  coupling  in  tetrahedron 
57,58,65,121,129  caused  a surface  sputtering  event. 

During  the  trace  of  the  sputtering 
mechanisms  at  600  eV,  the  third  and  fourth  layers  never 
contributed  an  active  atom  to  the  sputtering  mechanism; 
however,  these  layers  are  essential  to  the  mechanism  since 
the  coupling  of  momentum  in  the  Y direction  through  the 
tetrahedra  allows  the  squeeze  mechanism  to  occur  without 
driving  the  "squeezed"  atom  downward.  The  rising  of  rwo 
atoms  from  the  second  layer  was  responsible  for  seven  of  the 
eight  sputtering  events  observed,  and  one  event  was  caused 
by  the  lateral  movement  of  a second  layer  atom. 

2 • iC0j>eV  Results 

The  complete  41-impact- point  simulation  program  run 
with  bombarding  ion  energy  of  1 keV  required  183  minutes  of 
computer  time.  The  sputtering  ratio  obtained  with  the 
surface  binding  energy  set  at  1 eV  was  3.415,  with  140  atoms 
being  sputtered.  Figure  19  shows  the  atoms  sputtered  at 
this  energy.  The  higher  energy  bombardment  caused  an 
increase  in  the  number  of  atoms  sputtered  in  the  second  and 
third  layers  and  caused  a fourth-layer  atom  (atom  225)  to  be 
sputtered.  This  indicates  that  higher  energy  simulations 
must  not  only  increase  the  lattice  size  in  the  X and  2 
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directions,  bat  mast  add  more  crystal  layers  in  the  Y 
direction  in  order  to  adequately  contain  all  sputtering 
events.  Atom  225  was  sputtered  with  an  X velocity  component 
of  only  29  m/sec,  and  a Z velocity  component  of  86  mSec, 
while  the  normal  velocity  component  was  -3681  mSec. 
Investigation  of  this  type  event  may  show  that  the  atom  in 
fact  ransited  a (1G0)  channel  while  rising  to  the  surface. 
A normal  velocity  component  of  3861  m/sec  corresponds  to  a 
kinetic  energy  of  2.2  eV,  a low  energy  sputtering  event, 
which  may  be  below  the  final  value  determined  for  the 
surface  binding  energy  by  further  studies.  Figure  20  shows 
the  plot  of  the  1000  eV  sputtering  ratio  vs.  surface  binding 
energy.  Again  the  rapid  rise  of  the  sputtering  ratio  below 
3 eV  was  observed,  suggesting  a binding  energy  above  this 
value.  This  would  then  eliminate  many  of  the  lower  layer 
sputtering  events. 

3 • 6 CO-e V and  1 00 O^eV  Comparison 

Figure  21  shows  the  comparison  of  the  1000-eV  and 
600-eV  sputtering  ratios  plotted  vs.  surface  binding  energy. 
The  1000-eV  simulation  proved  to  have  a higher  sputtering 
ratio  at  each  value  of  3THR  as  might  be  expected,  but  the 
values  of  the  sputtering  ratios  closely  approached  each 
other  above  a binding  energy  of  8 eV.  Above  this  point  only 
the  highest  energy  sputtering  events  are  counted,  and  it  is 
reasonable  to  assume  the  value  of  ETHR  is  below  the  point 
where  the  sputtering  ratios  converge. 

3.  (110)  LATTICE  ORIENTATION 

The  (110)  lattice  orientation  proved  very  interesting 
because  the  sputtering  ratio  decreased  as  the  bombarding  ion 
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1 . oOOzil  Results 


The  600-eV  simulation  of  the  (110)  orientation 
utilized  139  minutes  of  computer  time  to  complete  30  impact 
points.  Shots  which  did  not  place  the  ion  near  the  channel 
center  required  approximately  the  same  number  of  timesteps 
as  the  (100)  orientation  computer  simulation,  while  shots  in 
which  the  bombarding  ion  transited  the  channel  typically 
required  fewer  timesteps.  The  600-eV  simulation  caused  60 
atoms  to  sputter  resulting  in  a sputtering  ratio  of  2.00  for 
a surface  binding  energy  of  1 eV.  Figure  23  shows  the 
sputtered  atoms  for  the  600-eV  simulation.  Again,  most 
sputtering  events  took  place  in  the  first  layer,  with  atom 
95  being  the  only  multiply  sputtered  atom  in  the  second 
la  yer . 


Figure  24 
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shows  the  plot  of  sputtering  ratio  vs. 
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normal  velocity  components  corresponding  to  an  energy  of  15 
eV  or  greater,  the  most  energetic  ejection  being  atom  95 
which  was  ejected  with  a normal  velocity  component  of 
-15,131  m/sec  when  atom  23  was  impacted  at  its  center  by  the 
bombarding  ion.  The  ratios  of  Y and  Z velocities  of  atom  95 
indicate  that  it  rose  1.81  LU  while  transiting  0.63  LU  in 
the  positive  Z direction.  Atom  95  then  rose  to  the  surface 
in  the  channel  formed  by  atoms  30,31,32,61,  and  63.  The 
(110)  channels,  therefore,  were  observed  to  allow  the  escape 
of  lower  atoms  with  high  energies,  while  focusing  their 
ejections  in  a near-normal  direction.  Atom  23,  the  next 
most  energetic  ejection  was  studied  in  the  mechanism  trace 
performed  at  the  most  active  impact  point. 

a.  Mechanis.  Summary 

The  most  active  impact  point  in  the  600  eV 
simulation  occured  where  the  ion  struck  atom  23  at  a point 
0.141  LU  to  the  right  of  center  and  0.2  LU  above  center  in 
the  positive  Z direction.  The  single-shot  simulation 
required  334  seconds  of  computer  time.  Atoms  4,14,18,23,43, 
and  61  were  sputtered.  Figure  25  illustrates  the  mechanisms 
described  in  the  individual  atom  traces. 

(1)  Atom  23.  The  Ar  ion  collided  with  atom  23 
at  timestep  1 and  imparted  negative  X and  Z velocity  to  it 
while  driving  it  downward.  At  timestep  20,  atom  23  reached 
its  point  of  closest  approach  to  atom  46,  passing  above  it 
at  a distance  of  0.54  LU.  A brief  interaction  with  atom  46 
slowed  a'tom  23*  s positive  Y velocity  slightly.  Atom  23 
continued  downward  and  in  the  negative  X and  Z direction, 
and  by  timestep  35  had  reached  its  closest  point  of  approach 
with  atom  45.  The  collision  with  atom  45  reversed  the 
descent  of  atom  23,  while  interactions  with  atoms  14  and  44 
prevented  direction  changes  in  the  Z direction.  At  timestep 
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60,  atom  23  was  no  longer  interacting  with  atoms  45,14,  or 
44,  but  was  interacting  only  with  atom  12.  The  collision 
with  atom  12  slowed  atom  23's  velocity  in  the  negative  X 
direction,  reversed  its  negative  Z velocity,  and  accelerated 
it  in  the  negative  Y direction.  At  timestep  80,  atom  23  was 
counted  as  a sputtered  atom.  At  that  tine,  atom  23  was 
still  interacting  with  atom  12,  but  the  computer  simulation 
ceased  accelerating  the  atom  because  it  had  already  risen 
0.4  LU  above  the  surface.  The  high  energy  ejection  of  atom 
23  was  a purely  surface  event,  with  the  impacted  atom 
traveling  a large  distance  across  the  surface  because  of  the 
vacancies  created  by  the  (110)  channels.  This  large  lateral 
displacement  is  not  possible  in  close  packed  orientations. 

(2)  Atom  lii*  Atom  23,  while  rebounding  from 
atom  45,  passed  within  0.618  LU  of  atom  14.  The  interaction 
of  atom  23  and  atom  14  caused  14  to  rise  with  a velocity  of 
-1013  m/sec.  At  timestep  50,  atom  14  was  squeezed  between 
atoms  23  and  15,  slowing  its  velocity  in  the  negative  Y 
direction.  Atom  14  then  continued  to  rise,  and  at  timestep 
60  was  interacting  with  no  other  atoms.  The  normal  velocity 
of  atom  14  at  that  time  was  -1421  n /sec,  insufficient  to 
cause  sputtering.  At  timestep  120,  atom  45,  which  was 
recoiling  upward  as  a result  of  its  previous  impact  by  atom 
23,  overtook  the  slowly  rising  atom  14.  The  resultant 
collision  slowed  atom  45,  but  increased  the  normal  component 
of  atom  14.  At  timestep  125  then,  atom  14  was  no  longer 
interacting  with  atom  45,  and  was  counted  as  a sputtered 
atom  with  a normal  velocity  component  of  -3600  m/sec.  The 
squeeze  mechanism  described  by  Harrison  et  al  [12]  was  the 
primary  mechanism  responsible  for  this  surface  sputtering 
event . 

(3)  Atom  4.  Atom  45,  after  being  struck  by 
atom  23,  moved  downward  under  atom  14.  At  timestep  70,  atom 
45  began  interaction  with  atom  36,  accelerating  36  upward 
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and  in  the  negative  X and  Z directions.  At  timestap  80, 
atom  45  reached  its  closest  point  of  approach  to  atom  36  and 
began  to  recoil  back  under  atom  14.  Atom  36  was  deflected 
slightly  by  atom  35,  and  by  timestep  110  had  moved  into 
interaction  with  atom  4.  Atom  4 was  accelerated  upward  by 
the  interaction  with  atom  36,  and  was  squeezed  slightly  by 
interaction  with  atom  3 at  timestep  115.  3y  timestep  125, 
atom  4 was  interacting  with  no  other  atoms,  and  was  rising 
with  a normal  velocity  component  of  -2674  m/sec.  At 
timestep  140,  atom  4 was  counted  as  a sputtered  atom.  The 
primary  mechanism  involved  was  the  mole  mechanism  described 
by  Harrison  et  al  [12].  This  is  distinguished  from  the 
scoop  or  squeeze  mechanism  by  the  low  angle  at  which  the 
approaching  atom  strikes  the  sputtered  atom. 

(4)  A tom  48.  Atom  1,  after  striking  atom  23, 
continued  into  the  second  layer.  Interaction  with  atoms 
113,119,  and  88  deflected  atom  1 into  the  (110)  channel  as 
it  passed  through  the  second  layer.  At  timestep  30,  atom  1 
reached  its  closest  point  of  approach  to  atom  83,  passing 
within  0.406  LTJ  and  accelerating  atom  88  downward  into 
collision  with  the  lower  part  of  atom  39.  Ey  timestep  45, 
atoms  88  and  89  reached  the  point  of  closest  approach  and 
atom  89  had  acquired  a Y velocity  component  of  -1309  m/sec. 
At  timestep  45,  atom  39  was  also  interacting  slightly  with 
atom  90,  and  was  being  deflected  in  the  positive  Z direction 
by  atom  48.  Atom  89  caused  atom  48  to  rise  from  the 
crystal.  By  timestep  80,  atom  48  was  no  longer  interacting 
with  any  other  atom,  and  at  timestep  145  it  was  counted  as  a 
sputtered  atom  with  a normal  velocity  component  of  -3400 
m/sec.  The  mechanism  involved  was  primarily  the  deep 
lechanKi. 

(5)  Atom  JJ.  After  being  deflected  by  atom  48, 
atom  89  collided  with  the  underside  of  atom  49.  Atom  49 
then  rose  in  the  lattice  pushing  atom  18  upward  from  the 
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surface.  At  timestep  100,  atoms  49  and  18  reached  the 
closest  point  of  approach,  their  separation  being  0.852  L3. 
At  timestep  105,  atom  18  had  acquired  a normal  velocity 
component  of  -3064  m/sec,  and  was  0.156  LU  above  the  lattice 
surface.  Atom  18  continued  to  rise  without  any  other 
interactions,  and  was  counted  as  a sputtered  atom  at 
timestep  125.  The  primary  mechanism  was  again  the  deep 
mechanism  caused  by  the  rising  of  atom  89  from  the  second 
layer . 


(6)  Atom  §2*  Atom  1,  when  transiting  the 
second  layer,  caused  the  outward  movement  of  atoms  118,119, 
and  83.  Atom  118  moved  under  atom  95, and  at  timestep  70, 
collided  with  atom  94.  Atom  94  moved  upward,  and  while 
still  interacting  with  atom  118,  caused  atom  61  to  rise  from 
the  surface.  At  timestep  90,  atom  94  had  reached  a maximum 
upward  velocity  of  -3106  m/sec,  and  atom  61  had  reached  an 
upward  velocity  of  -1663  m/sec.  By  timestep  110,  the 
coupling  of  momentum  from  atom  118  through  94  to  atom  61  was 
complete,  with  atom  61  rising  at  -4071  m/sec,  and  atom  94 
having  slowed  to  -731  m/sec.  Atom  61  continued  to  rise  with 
no  further  interactions,  and  was  counted  as  a sputtered  atom 
at  timestep  145.  The  mechanism  was  the  deep  mechanism 
caused  by  atom  94. 


Six  sputtering  events  occured  at  this 
impact  point;  three  sputtering  events  were  the  result  of  the 
deep  mechanism,  two  sputtering  events  were  surface  events 
characterized  by  the  squeeze  or  mole  mechanism,  and  the  most 
energetic  event  was  characterized  by  a mechanism  not  seen  in 
the  more  close-packed  lattice  orientations.  This  large 
lateral  displacement  followed  by  one  or  more  surface 
reflections  can  perhaps  best  be  described  as  a skip 
mechanism.  The  other  two  surface  mechanisms  observed  (the 
sputtering  of  atoms  4 and  14)  were  the  indirect  result  of 
the  "skipping"  of  atom  23. 
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2. 


1 GO Q- eV  Results 

The  30-impact-point,  1-keV  simulation  program 
required  141  minutes  of  computer  time,  two  minutes  longer 
than  the  600  eV  simulation.  With  the  binding  energy  set  at 
1 eV,  56  atoms  were  sputtered,  resulting  in  a sputtering 
ratio  of  1.866.  The  lower  sputtering  ratio  was  unexpected 
at  a bombarding  energy  of  1 keV ; therefore,  the  mechanisms 
of  sputtering  were  again  traced  at  the  same  impact  point 
that  was  examined  at  600  eV.  Figure  26  shows  the  atoms 
sputtered  in  the  1-keV  simulation.  At  1 keV,  the  sputtering 
was  again  primarily  a surface  phenomenon,  with  only  three 
single  sputs  occuring  in  the  second  layer,  and  no  atoms 
being  sputtered  from  any  layer  lower  than  the  second.  The 
1-keV  simulation  caused  the  sputtering  of  five  boundary 
atoms,  while  the  600-eV  simulation  sputtered  four  atoms  on 
the  surface  layer  boundary.  Only  two  atoms  were  ejected 
with  normal  velocity  components  corresponding  to  kinetic 
energies  greater  than  15  eV,  as  compared  with  4 high  energy 
ejections  in  the  600-eV  simulation.  At  impact  point  0 (dead 
center  of  atom  23) , atom  95  was  ejected  with  a normal 
velocity  component  of  -16,296  m/sec. ^ Atom  95  was  ejected 
through  the  same  channel  as  it  was  in  the  600-eV  simulation. 
Figure  27  shows  the  plot  of  sputtering  ratio  vs.  ETHR  for 
the  1 keV  simulation.  Again,  a rapid  rise  in  the  sputtering 
ratio  was  observed  below  3 eV.  The  decrease  between  4 eV 
and  8 eV,  however,  was  not  as  gently  sloping  as  the  600-eV 
run.  The  sputtering  ratio  continued  to  drop  fairly  rapidly 
until  9 eV,  at  which  point  only  the  two  highest-en ergy  atoms 
were  sputtered. 
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a.  Mechanism  Summary 

The  1-keV  single-impact- point  program  required 
339  seconds  of  computer  time,  5 seconds  longer  than  the  600 
eV  simulation.  Atoms  14,18,23,48,  and  65  were  sputtered,  a 
decrease  of  1 atom  from  the  600  eV  simulation.  Two  atoms 
(atoms  4 and  61)  which  were  sputtered  at  600  eV,  were  not 
sputtered  at  1 keV.  Atom  65  was  sputtered  at  1 keV,  but  not 
in  the  lower-energy  simulation.  The  mechanism  trace  was 
performed  and  compared  at  each  step  with  the  corresponding 
600-eV  program.  Figure  28  illustrates  the  individual  atom 
traces . 

(1)  Atom  23.  At  timestep  1,  atom  23  was 
impacted  0.141  LU  to  the  right  of  center  and  0.2  LU  above 
center  in  the  Z direction.  The  impact  drove  atom  23 
downward,  and  in  the  negative  X and  Z directions.  Atom  23 
was  deflected  slightly  by  atom  46,  and  at  timestep  45  had 
been  deflected  upward  by  atom  45.  At  timestep  70,  atom  23 
was  interacting  with  atoms  12  and  13,  striking  slightly 
below  these  atoms.  At  timestep  80,  atom  23  reached  the 
closest  point  of  approach  to  to  atom  12,  and  was  deflected 
upward  in  a nearly  normal  direction.  Atom  12  continued 
upward,  interacting  primarily  with  atom  12,  and  at  timestep 
100  had  passed  upward  out  of  interaction  range  with  any 
atom.  Atom  23  continued  upward  with  a normal  velocity 
component  of  -7500  m/sec  and  nearly  zero  X and  Z velocities 


and  was  counted  as  a sputtered  atom  at  timestep  120.  The 
skip  mechanism  again  caused  the  sputtering  of  atom  23,  as  it 
did  in  the  600  eV  simulation;  however,  atom  23  was  not 
deflected  upward  as  much  by  atom  45  and  therefore  impacted 
slightly  below  atom  12  and  13,  slowing  its  normal  velocity 
component  rather  than  accelerating  it  as  was  the  case  in  the 
600  eV  simulation. 
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(2)  Atom  14.  Atom  23,  while  rebounding  from 
atom  45,  passed  within  0.601  LU  of  atom  14.  The  interaction 
of  atom  14  and  23,  combined  with  slight  recoils  from  atoms 
13  and  15,  caused  atom  14  to  rise  from  the  lattice.  Atom  14 
continued  to  rise,  and  at  timestep  105  was  interacting  with 
no  other  atoms.  At  timestep  115,  atom  14  was  counted  as  a 
sputtered  atom.  Atom  14  possessed  a normal  velocity 
component  of  -4491  m/sec  at  that  time.  Once  again  the  scoop 
or  squeeze  mechanism  was  responsible  for  a surface 
sputtering  event.  In  this  case  14  was  not  impacted  by  atom 
45  after  being  squeezed  by  atom  23,  but  received  sufficient 
velocity  from  the  primary  mechanism  itself.  After  being 
driven  downward  by  atom  23,  atom  45  came  to  rest  0.711  LU 
below  its  initial  position. 

(2)  Atom  43.  Atom  1,  after  striking  atom  23, 
continued  into  the  second  layer.  At  timestep  30,  atom  1 
impacted  atom  88,  causing  atom  88  to  move  in  the  positive  X 
direction  underneath  atom  89.  At  timestep  45,  atom  89 
started  to  rise  from  the  second  layer,  and  by  timestep  50 
atom  89  was  accelerating  atom  48  upward.  At  timestep  70, 
atom  89  reached  its  closest  point  of  approach  to  atom  43, 
rising  to  within  0.765  LU.  Atom  89  was  slowed  by  atom  48, 
and  by  timestep  95  was  traveling  downward  toward  its  initial 
position.  Atom  48  continued  to  rise,  interacting  slightly 
with  atom  47,  and  pushing  atom  25  upward  at  a velocity  too 
low  to  cause  it  to  sputter.  By  timestep  95,  atom  48  was 
interacting  with  no  other  atom  and  was  rising  with  a normal 
velocity  component  of  -3552  m/sec.  At  timestep  155,  atom  U8 
had  risen  0.418  LU  above  the  surface  and  was  counted  as  a 
sputtered  atom.  The  mechanism  involved  was  essentially  the 
sane  deep  mechanism  involved  in  the  sputtering  of  atom  48 
during  the  600-eV  simulation.  The  only  difference  was  that 
atom  88  was  impacted  more  directly  by  atom  1,  because  the 
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(4)  Atom  13..  At  timestep  90,  atom  89  after 
being  deflected  by  its  collision  with  atom  48,  interacted 
with  atom  49.  At  timestep  100,  atom  49  having  been 
accelerated  in  the  positive  X and  negative  Z directions, 
rose  under  atom  13.  At  timestep  120,  atom  49  was  at  its 
closest  point  of  approach  to  atom  13,  passing  under  atom  18 
at  a distance  of  0.836  LU.  Atom  18  rose  without  interacting 
with  any  other  atom,  and  at  timestep  125  possessed  zero 
potential  energy.  At  timestep  135,  atom  18  was  counted  as  a 
sputtered  atom,  possessing  a normal  velocity  component  of 
-4230  m/sec.  The  primary  mechanism  involved  was  the  deep 
mechanism  caused  by  the  rising  of  atom  89  from  the  second 
layer,  just  as  was  the  case  in  the  600-eV  sputtering  event. 


(5)  Atom  65.  Atom  1,  when  colliding  with  atom 
88  in  second  layer,  also  interacted  with  atom  110, 
accelerating  it  in  the  positive  X and  Z direction  and 
driving  it  downward.  At  timestep  45,  atom  119  moved  within 
interaction  range  with  atom  120,  accelerating  atom  120  in 
the  positive  X and  Z direction  with  a virtually  negligible  Y 
velocity  component.  At  timestep  65,  atom  120  interacted  not 
only  with  atom  119,  but  experienced  grazing  collisions  with 
atoms  97,121,  and  161  (atom  161  is  located  directly  below 
atom  97  in  the  fourth  layer) . Atom  120  was  focused  then 
into  collision  with  atom  98,  and  was  within  interaction 
range  at  timestep  100.  At  timestep  100,  atom  119  had 
recoiled  slightly,  and  was  no  longer  interacting  with  atom 
120.  Atom  120  was  rising  slightly  as  a result  of  its 
collision  with  atom  161,  and  rose  underneath  atom  98.  Atom 
98  then  rose  toward  the  first  layer,  and  at  timestep  115 
began  interacting  with  atom  65.  Atom  65  rose  from  the 
surface  interacting  with  only  atom  98.  At  timestep  130, 
atom  65  had  acquired  a normal  velocity  component  of  -3246 
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m/sec,  and  was  interacting  with  no  other  atoms.  At  timestep 
130,  atom  65  reached  a height  of  0.499  LU  above  the  surface 
and  was  counted  as  a sputtered  atom.  Atom  65  was  sputtered 
by  the  deep  mechanism.  Essentially  the  same  mechanism 
occured  in  the  600  eV  simulation,  but  the  lower  energy  Ar 
ion  was  deflected  into  the  (110)  channel  and  therefore  did 
not  transfer  sufficient  momentum  to  atom  119.  Atom  98  did, 
in  fact,  rise  at  a much  lower  velocity  causing  atom  65  to 
rise.  At  program  termination,  however,  atom  65  was  rising 
with  a normal  velocity  component  of  -877  m/sec  and  was  still 
0.2  LU  below  the  lattice  surface. 

Atoms  61  and  4 were  sputtered  in  the  600-eV 
simulation,  but  not  in  the  1-keV  simulation.  Atom  4 was  not 
sputtered  because  atom  45  was  driven  deeper  into  the 
crystal,  and  came  to  rest  in  a lower  layer,  rather  than 
recoiling  into  atom  13  as  it  did  in  the  600-eV  simulation 
(Fig.  25).  Atom  61  was  not  sputtered  because  the  Ar  ion  was 
not  deflected  into  the  (110)  channel.  Atom  118,  therefore, 
did  not  achieve  a high  outward  velocity.  In  the  1-keV 
simulation,  atom  118  moved  into  the  vacancy  left  dv  atom  119 

i 

rather  than  in  the  negative  X direction.  Atom  94  never 
possessed  sufficient  energy  to  be  listed. 

3 • 2 CO Q-eV  Results 

The  2-keV  30- impact-point  program  required  135 
minutes  of  computer  time,  4 minutes  less  than  the  600  aV 
simulation.  With  ETHR  set  at  1 eV,  44  atoms  were  sputtered, 
resulting  in  a sputtering  ratio  of  1.467.  This  sputtering 
ratio  was  lower  than  the  sputtering  ratio  obtained  in  the 
1-keV  simulation9  Figure  29  shows  the  sputtered  atoms 
obtained  at  1 eV  binding  energy.  Once  again,  the  sputtering 
was  observed  to  be  predominantly  a surface  phenomenon,  with 
two  sputtering  events  occuring  in  the  second  layer  (atoms  95 
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and  127),  and  a single  spattering  event  occuring  in  the 
fourth  layer  (atom  223)  . Two  high  energy  ejections  were 
observed,  occuring  at  the  same  impact  points  as  the  high 
energy  ejections  in  the  1-keV  program.  Atom  95  was  ejected 
with  a normal  velocity  component  of  -16,020  m/sec,  176  m/sec 
slower  than  the  1 keV  event,  and  atom  63  was  ejected  at 
-17,835  m/sec,  2,353  m/sec  more  than  the  1-keV  event.  Atom 
95  again  exited  the  surface  through  the  (110)  channel,  as  it 
did  in  the  two  lower-energy  simulations.  Atom  223,  the 
fourth-layer  sputtering  event  also  exited  the  surface 
through  the  same  (110)  channel.  Only  two  boundary  atoms 
were  sputtered  in  the  surface  layer,  as  compared  to  five  in 
the  1-keV  simulation,  and  four  in  the  600-eV  simulations; 
however,  the  sputtering  of  a fourth-layer  atom  showed  that 
not  only  must  the  crystal  be  expanded  in  the  X and  Z 
directions,  but  more  layers  must  be  added  if  containment  of 
the  events  is  to  be  achieved.  The  fourth-layer  sput  occured 
at  an  energy  of  2.38  eV,  again  a relatively  low  energy 
sputtering  event.  Figure  30  shows  the  plot  of  sputtering 
ratio  vs.  surface  binding  energy.  The  rapid  rise  of 
sputtering  ratio  below  3 eV  is  again  evident,  and  a plateau 
value  of  0.133  is  reached  at  8 and  9 eV,  similar  to  the 
1-keV  case.  A surface  oinding  energy  above  the  3-eV  point 
would  eliminate  the  lower-layer  sputtering  event. 


a.  Mechanism  Summary 


The  2-keV  single-shot  simulation  required  369 
seconds  of  computer  time,  30  seconds  longer  that  the  1 keV 
simulation.  Five  first- layer  atoms  were  sputtered,  the  same 
number  as  in  the  1 keV  simulation;  however,  atoms  23,43,  and 


65  were  not 

sputtered  at  2 keV, 

and  atoms  4, 

12, 

and  13 
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sputtered 

instead,  indicating 

another  change 

of 

the 

predominant 

mechanisms.  The 
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in  figure  31.  The 
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sputtered  by  the  skip  mechanism  in  both  the  600 -eV  and  1-keV 
simulations,  but  not  sputtered  in  the  2-keV  simulation  is 
also  included. 

(1)  Atom  23.  At  timestep  1,  atom  23  was 
impacted  by  the  Ar  ion.  Atom  23  moved  over  atom  46  and 
collided  with  atom  45  moving  atom  45  downward  into  the 
crystal.  Atom  23  then  rebounded  from  atom  45,  and  struck 
atoms  12  and  13  much  as  it  did  in  the  600-eV  and  1-keV 
programs;  however,  atom  23  struck  below  these  two  atoms  and 
was  deflected  downward.  Atom  23  was  deflected  further 
downward  by  collision  with  atom  35,  and  continued  down  into 
the  second  layer.  At  timestep  160,  atom  23  exited  the 
crystal  at  its  left  side  while  traveling  between  the  second 
and  third  layers,  carrying  with  it  21.29  eV  of  kinetic 
energy  after  having  lost  '116.44  e7  of  its  initial  kinetic 
energy  by  collisions  within  the  lattice.  The  skip  mechanism 
seems,  then,  to  be  energy  dependent.  At  higher  energies  the 
"skipping"  atom  causes  too  much  surface  damage  when 
colliding  with  surface  atoms,  therefore  moving  into  the 
crystal  rather  than  rebounding  across  its  surface.  The 
absence  of  this  mechanism  at  higher  energies  could,  in  part 
at  least,  account  for  the  decreased  sputtering  ratio  at 
higher  energies. 

(2)  Atom  1j*.  At  timestep  50,  atom  23,  while 
still  traveling  downward  over  atom  45,  began  interacting 
with  atom  14.  Atom  14  was  forced  upward  from  the  lattice  by 
this  interaction.  At  timestep  75,  atom  23  reached  its  point 
of  closest  approach  to  atom  14,  passing  below  it  at  a 
distance  of  0.513  LO.  At  the  same  time  atom  14  experienced 
grazing  collisions  with  atoms  13  and  15,  deflecting  it 
upward  even  further.  At  timestep  110,  atom  14  had  risen 
0.40  LO  above  the  surface.  Although  it  was  still 
interacting  with  atom  23,  it  was  counted  as  a sputtered 
atom.  The  mechanism  which  caused  atom  14  to  sputter  was  no 
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longer  predominantly  the  squeeze  or  scoop  mechanism,  but 
could  be  classified  as  the  mole  mechanism  because  of  the 
lower  impact  on  atom  14. 

(3)  Atom  12.  At  timestep  75,  atom  14 
interacted  with  atom  13  causing  it  to  move  -0.007  L(J  in  the 
X direction,  and  -0.006  LU  in  the  Z direction,  while  rising 
-0.004  LU.  Atom  13  continued  to  move  in  the  negative  X,Y, 
and  Z directions  at  a very  low  velocity  until  atom  23  began 
interaction  with  it.  At  timestep  100  atom  13  was  passed  by 
atom  23  at  a distance  of  0.621  LU.  At  the  closest  point  of 
approach,  atom  23  was  0.171  LU  below  atom  13,  thereby 
accelerating  atom  13  upward  from  the  surface.  Atom  13 
continued  to  rise,  interacting  with  no  other  atoms,  and  at 
timestep  130  was  counted  as  a sputtered  atom.  The  low 
impact  of  atom  23  on  atom  13  classifies  this  as  a mole 
mechanism  also. 

(4 ) Atom  \2.  After  passing  atoms  14  and  13, 
atom  23  collided  with  the  underside  of  atom  12  and  was 
sharply  deflected  downward.  Atom  23  moved  into  interaction 
range  with  atom  12  at  timestep  100,  and  reached  its  closest 
point  of  approach  at  timestep  105,  passing  under  atom  12  at 
a distance  of  0.508  LU.  Atom  12  was  accelerated  in  the  X 
direction  tc  a velocity  of  -16,863  m/sec,  and  upward  with  a 
velocity  of  -10,307  m/sec.  At  timestep  115,  atom  12  was  no 
longer  interacting  with  any  atom,  and  it  was  counted  as  a 
sputtered  atom.  Once  again  a surface  atom  was  sputtered  by 
the  mole  mechanism. 

(5)  Atom  4.  After  being  deflected  downward  by 
atom  12,  atom  23  continued  in  the  negative  X direction  and 
at  timestep  120  moved  into  interaction  range  with  atom  35. 
At  timestep  125,  atom  23  reached  its  closest  point  of 
approach  to  atom.  35,  accelerating  it  upward  into  collision 
with  atom  4.  Atom  4 rose  from  the  lattice  interacting  with 
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no  other  atom,  and  at  timestep  140  had  moved  out  of 
interaction  range  with  atom  35.  At  timestep  150,  atom  4 was 
counted  as  a sputtered  atom,  possessing  a normal  velocity 
component  of  -3,111  m/sec.  The  mole  mechanism  caused  the 
transfer  of  momentum  from  the  center  atom  of  the  imcomplete 
tetrahedron  4,35,36,  and  68  to  the  upper  atom  (atom  4), 
causing  another  surface  sputtering  event.  Atom  4 was  not 
sputtered  in  the  1-keV  simulation,  but  was  sputtered  in  the 
600-eV  simulation.  In  the  600  eV  simulation,  however,  the 
recoil  of  atom  45  caused  atom  36  to  impact  atom  4 and  raise 
it  from  the  surface.  In  the  2-keV  case,  a totally  different 
interaction  chain  was  involved. 

(6)  Atom  1 6 . After  striking  atom  23,  atom  1 
moved  into  the  second  layer  striking  atom  88.  Atom  88 
traveled  in  the  positive  X direction,  and  rose  from  its 
position  in  the  second  layer,  forcing  atom  89  downward  and 
in  the  positive  X direction.  At  timestep  80,  atom  89 
interacted  with  atom  48,  causing  it  to  rise  slowly  from  the 
crystal  surface.  Atom  89  also  recoiled  from  atom  112,  and 
rose,  interacting  with  atom  49  at  timestep  130.  At  timestep 
135,  atom  18  was  lifted  from  the  crystal  surface  by  the 
rising  atom  49,  and  interacted  briefly  with  atom  17,  which 
was  moving  as  the  result  of  a collision  with  atom  15.  At 
timestep  140,  atom  18  was  interacting  only  with  atom  49.  At 
timestep  145,  atom  13  was  interacting  with  no  other  atoms 
and  was  rising  from  the  crystal  with  a normal  velocity 
component  corresponding  to  a kinetic  energy  of  1.01  eV.  At 
timestep  160,  atom  18  was  counted  as  a sputtered  atom.  The 
mechanism  was  again  the  deep  mechanism  caused  by  the  rising 
of  atom  89  from  the  second  layer. 

In  both  the  600-eV  and  1-keV  simulations, 
atom  48  was  sputtered  by  the  deep  mechanism  caused  by  the 
rising  of  atom  89.  In  the  2 keV  simulation,  however,  atom 
88  rose  above  atom  89  rather  than  squeezing  it  upward,  and 
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atom  89  rose  as  the  result  of  recoil  from  atom  112.  This 
lower  energy  rebound  did  not  accelerate  atom  48  sufficiently 
to  cause  it  to  be  counted  as  a sputtered  atom.  In  the  1-keV 
simulation,  atom  65  was  sputtered  by  atom  98  as  it  rose  from 
the  second  layer.  In  the  2-keV  simulation,  atom  119  was  not 
impacted  by  atom  1,  and  remained  essentially  at  rest  in  the 
crystal;  therefore,  atom  120  did  not  transfer  a significant 
amount  of  energy  tc  atom  93. 

Although  atom  65  did  rise  slightly,  it  was  still  below  the 
crystal  surface  at  program  termination. 

4 . Comparison 


A comparison  of  the  sputtering  ratios  vs.  ET HR  for 
each  bombarding  energy  yielded  no  insight  into  the  value  for 
surface  binding  energy  of  the  (110)  orientation  of  silicon. 
Figure  31  shows  that  for  a binding  energy  below  3-eV,  the 
sputtering  ratios  of  the  600-eV  program  were  higher  than  the 
1-keV  or  2-keV  programs.  Between  3 eV  and  7 eV,  the  1-keV 
program  possessed  the  highest  sputtering  ratios,  and  between 
8 eV  and  15  eV,  the  600  eV  sputtering  ratios  were  again  the 
highest.  At  a binding  energy  of  7 eV,  all  three  programs 
displayed  the  same  sputtering  ratio.  Comparison  of  the 
mechanism  traces  showed  that  the  skip  mechanism  probably 
occurs  only  at  low  bombarding  energies,  and  also  that  as 
bombarding  energy  increases,  surface  damage  caused  by 
driving  surface  atoms  deeper  into  the  crystal  decreases  the 
number  of  surface  sputtering  events. 

C.  (Ill)  LATTICE  ORIENTATION 

The  (111)  lattice  orientation,  being  a more  open 
orientation  than  the  (100),  was  expected  to  have  a lower 
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sputtering  ratio  at  corresponding  energies.  This  was 
confirmed  at  both  600  eV  and  1 keV.  The  sputtering  ratio  at 
1 keV  was  higher  than  the  sputtering  ratio  at  600  e'/, 
showing  that  the  unexpected  results  encountered  in  the  (110) 
orientation  were  not  duplicated  in  the  (111)  orientation. 
Figure  33  shows  isometric  projection  of  the  (111) 
orientation.  The  fourth  layer  is  a duplicate  of  the  first, 
and  was  eliminated  from  the  figure  in  the  interest  of 
legibility. 


1 . cOO-eV  Results 


The  complete  36-impact-point  600-eV  simulation 
required  139  minutes  of  computer  time.  With  the  binding 
energy  set  at  1 eV , 51  atoms  were  sputtered  yielding  a 
sputtering  ratio  of  1.417.  All  sputtered  atoms  were  located 
in  the  first  three  layers,  with  four  multiple  and  three 
single  sputs  originating  in  the  second  layer,  and  three 
single  sputtering  events  originiating  in  the  third  layer 
(Fig.  34).  Nine  atoms  located  on  the  X and  Z boundaries 
were  sputtered,  indicating  that  the  lattice  dimensions  were 
too  small  to  sufficiently  contain  all  events  at  600  eV. 
Four  atoms  were  ejected  at  energies  greater  than  15  eV. 
Three  of  the  high  energy  ejections  were  surface  layer  events 
(atoms  24,27,59),  and  one  ejection  (atom  171)  originated  in 
the  third  layer.  Atom  171  was  sputtered  from  the  third 
layer  with  a normal  velocity  component  corresponding  to  28 
eV,  indicating  that  the  open  lattice  allows  channeling  of 
these  third  layer  atoms  to  the  surface  without  experiencing 
hard  collisions  enroute. 

Figure  35  shows  the  plot  of  the  600-eV  sputtering 
ratio  vs.  ETHR.  Again,  the  rapid  rise  of  the  sputtering 
ratio  below  3 eV  was  observed,  suggesting  that  the  bindinq 
energy  is  above  3 eV. 
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Mechanism  Summary 


a . 


A five-minute,  sin gle- im pact- point  program  was 
run  with  600-eV  bombarding  energy,  and  the  impact  point 
displaced  0.141  LU  from  the  impact  atom  center  in  the  X 
direction  and  0.133  LU  from  the  atom  center  in  the  Z 
direction.  At  this  active  impact  point,  six  atoms  were 
sputtered  from  the  surface  layer,  and  two  atoms  were 
sputtered  from  the  second  layer.  Figure  36  illustrates 
traces  of  the  atoms  listed  separately  in  the  summary. 

(1)  Atom  30.  At  timestep  1,  atom  39  was 
impacted  by  the  Ar  ion.  Atom  39  moved  downward  while 
traveling  in  the  negative  X and  Z directions,  and  at 
timestep  25  moved  into  interaction  range  with  atom  30.  Atom 
30  was  accelerated  downward  at  13,328  m/sec,  and  at  timestep 
35,  collided  with  atom  85.  Atom  30  recoiled  upward  from 
atom  35  and  rose  beneath  atom  24.  At  timestep  65,  atom  30 
passed  within  0.651  LU  of  atom  24  and  was  deflected  in  the 
positive  Z direction.  After  passing  beneath  atom  24,  atom 
30  collided  with  atom  32  and  was  deflected  upward.  At 
timestep  95,  the  rapidly  rising  atom  30  interacted  with  atom 
31,  causing  its  normal  velocity  component  to  be  slowed  from 
-5729  m/sec  to  -4405  m/sec.  At  timestep  110,  a brief 
interaction  with  atom  23  further  slowed  atom  30's  normal 
velocity  component  to  -2845  m/sec.  At  timestep  125,  atom  30 
had  risen  0.119  LU  above  the  lattice  surface,  and  was 
interacting  with  no  other  atoms.  At  timestep  140,  atom  30 
was  counted  as  a sputtered  atom.  This  interaction  was  the 
result  of  a second-layer  reflect  ion, the  only  pure  reflection 
observed  in  the  mechanism  traces  performed. 
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The  placement  of  the  atoms  in  the  second  and  third  layers 
with  "pseudo-channels'’  to  the  surface  makes  this  type  of 
event  more  probable  in  the  (111)  orientation  than  in  the 
(100)  or  (110)  orientations. 

(2)  Abo™  24.  Atom  30,  while  rising  to  the 
surface,  began  interaction  with  atom  24  at  timestep  55. 
Atom  30  passed  below  atom  24  at  a distance  of  0.651  L'J, 
accelerating  it  in  the  negative  Y direction.  Atom  24  was 
also  accelerated  in  the  positive  X and  negative  Z directions 
by  atom  30,  and  at  timestep  95  it  was  slowed  by  collision 
with  atom  16.  At  timestep  105,  atom  24  interacted  with  both 
atoms  15  and  16,  and  was  deflected  toward  atoms  25  and  25. 
The  interaction  with  atoms  25  and  26  at  timestep  120  slowed 
atom  24 ' s normal  velocity  component  to  -2071  m/sec,  reducing 
its  energy  below  the  1 eV  required  to  sputter.  At  timestep 
130,  atom  24  was  still  0.173  LU  below  the  surface,  and  its 
normal  velocity  component  had  fallen  to  -1255  m/sec.  At 
timestep  155,  atom  79  rose  under  atom  24  as  the  result  of 
impact  by  atom  85.  Atom  74  accelerated  atom  24  upward,  and 
by  timestep  160  atom  24  was  interacting  with  no  other  atom 
and  was  rising  with  sufficient  energy  to  be  counted  as  a 
sputtered  atom.  At  timestep  165,  atom  24  was  counted  as  a 
sputtered  atom.  The  primary  mechanism  which  caused  atom  24 
to  sputter  was  the  deep  mechanism  caused  by  atom  79  rising 
from  the  second  layer.  This  was  assisted  by  the  reflection 
of  atom  30  from  the  second  layer. 
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(4)  Atom  25.  At  timestep  120,  atom  24  was 
slowed  by  atom  25.  The  interaction  with  atom  24  caused  atom 
25  to  rise  from  the  crystal  surface.  At  timestep  130,  atom 
25  had  risen  0.292  LU  above  the  lattice  surface  and  was  no 
longer  interacting  with  atom  24.  At  timestep  135,  atom  25 
was  counted  as  a sputtered  atom.  Atom  24  struck  atom  25  at 
a point  which  distinguishes  this  event  as  a mole  mechanism. 

(5)  Atom  47.  The  Ar  ion  (atom  1)  passed 
downward  into  the  second  layer  after  colliding  with  atom  39. 
During  its  passage  downward,  it  passed  within  0.671  LU  of 
atom  48,  causing  atom  48  to  rise  upward  under  atom  47.  At 
timestep  35,  atom  48  caused  atom  47  to  rise  upward.  Atom  47 
was  accelerated  to  a normal  velocity  of  -2880  m/sec  by  atom 
48,  and  rose  from  the  surface  interacting  with  no  other 
atom.  At  timestep  95,  atom  47  was  counted  as  a sputtered 
atom.  This  type  of  mechanism  shows  again  the  interaction  of 
nearest  neighbors  in  the  basic  tetrahedra,  and  can  best  be 
characterized  by  the  deep  mechanism. 

(6)  Atom  119.  The  Ar  ion,  after  passing 
through  the  second  layer,  passed  within  0.391  LU  of  atom 
175,  accelerating  it  in  the  positive  Z direction.  The 
movement  of  atom  175  caused  atom  120  to  rise  under  atom  119. 
At  timestep  60,  atom  119  was  accelerated  upward.  At 
timestep  70,  atom  119  was  deflected  slightly  by  atom  58,  and 
at  timestep  75  it  was  deflected  slightly  by  atom  56.  At 
timestep  85,  atom  119  was  0.935  LU  below  the  surface  and  was 
rising  freely,  interacting  with  no  other  atom,  with  a normal 
velocity  component  of  -3861  m/sec.  At  timestep  110,  atom 
119  passed  within  0.797  LU  of  atom  48,  and  was  slowed  to 
-3747  m/sec.  By  timestep  120,  atom  119  was  0.384  LU  below 
the  surface,  and  was  again  rising  freely  with  a normal 
velocity  component  of  -3645  m/sec.  At  timestep  145,  atom 
119  was  counted  as  a sputtered  atom.  The  primary  sputtering 
mechanism  was  the  mole  interaction  of  atoms  175  and  120. 
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(7)  Atom  63.  The  interaction  of  atoms  120  and 
64  caused  atom  64  to  rise  below  atom  63.  At  timestep  80, 
atom  64  pushed  atom  63  upward  from  the  surface.  Atom  63 
interacted  with  no  other  atoms,  and  at  timestep  105  was 
counted  as  a sputtered  atom.  The  mechanism  was  the  deep 
mechanism  caused  by  atom  120. 

(8)  Atom  95.  After  colliding  with  atom  39, 
atom  1 continued  downward  into  the  second  layer.  At 
timestep  25,  atom  1 passed  within  0.404  LO  of  atom  103, 
accelerating  atom  103  into  collision  with  atom  98.  Atom  103 
passed  within  0.508  LU  of  atom  98,  accelerating  atom  98 
downward  and  in  the  negative  X and  2 directions.  At 
timestep  115,  atom  98  entered  into  interaction  with  atom  95, 
accelerating  atom  95  upward,  and  in  the  negative  X and  Z 
directions.  at  timestep  120,  atom  98  passed  its  closest 
point  of  approach  to  atom  95.  At  timestep  130,  atom  95  had 
risen  -0.175  LU,  and  was  no  longer  interacting  with  any 
atom.  Atom  95  continued  to  rise,  and  at  timestep  175  was 
counted  as  a sputtered  atom.  The  primary  sputtering 
mechanism  was  the  mole  mechanism  which  occurred  when  atom  98 
collided  with  the  lower  part  of  atom  95. 


The  mechanisms  observed  in  the  (111) 
orientation  were  all  of  the  type  observed  by  Harrison  et  al 
[12].  Harrison  et  al  cited  rare  examples  of  the  reflection 
of  impacted  atoms  from  lower  layers  resulting  in  a 
sputtering  event.  Tne  open  nature  of  the  (111)  orientation 


makes  this  type 


event  more  probable  in  the  diamond 


lattice  than  in  the  close-packed  fee  crystal  investigated  by 
Harrison  et  al. 
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2 • 1 00  O-eV  Results 

The  1 keV  simulation  program  required  162  minutes  of 
computer  time,  23  minutes  longer  than  the  600-eV  simulation. 
A total  of  56  atoms  were  sputtered  with  the  binding  energy 
set  at  1 keV,  yielding  a sputtering  ratio  of  1.556.  Figure 
37  shows  the  atoms  sputtered  at  this  energy.  All  sputtering 
events  occured  in  the  first  three  crystal  layers;  however, 
ten  atoms  were  spattered  in  the  second  layer  as  compared 
with  seven  in  the  600-eV  simulation,  and  six  third-layer 
atoms  were  sputtered,  twice  as  many  as  in  the  600-eV 
simulation.  Eleven  boundary  atoms  were  sputtered,  one  more 
than  in  the  600  eV  simulation.  Four  atoms  were  ejected  with 
normal  velocity  components  corresponding  to  a kinetic  energy 
greater  than  15  eV . Three  of  the  high-energy  sputs  were 
surface  events,  while  one  high-energy  sput  originated  in  the 
third  lattice  layer.  The  third-layer  sput  (atom  179) 
traveled  in  the  negative  X direction  while  rising  to  the 
surface.  Since  atom  179  was  located  at  the  left  edge  of  the 
crystal,  it  traveled  outside  the  lattice  boundary  while 
rising,  and  might  have  been  slowed  if  the  lattice  were 
larger . 

Figure  38  shows  the  plot  of  sputtering  ratio  vs. 
ETHR  for  the  1 keV  simulation.  The  familiar  rise  below  3 eV 
binding  energy  was  again  evident.  The  rate  of  decrease  of 
sputtering  ratio  above  3 eV  was  approximately  the  same  as 
the  600  eV  simulation's  rate  of  decrease.  The  sputtering 
ratio  reached  a plateau  of  0.250  at  9 and  10  eV,  and 
decreased  to  0.111  at  15  eV. 

The  increased  sputtering  of  deeper  atoms  as  the 
bombarding  energy  increased  was  expected  because  of  the 
geometry  of  the  (111)  orientation.  The  open  appearance  of 
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the  "pseudo-channels"  is  terminated  in 
rising  atoms  below  this  depth  enc 
crystal.  Considerable  damage  to  the  s 
would  be  required  to  create  vacancies 
events. 
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value  of  binding  energy  is  probably  below  this  value.  The 
number  of  boundary  atoms  sputtered  at  600  eV  and  at  1 keV 
indicated  that  the  lattice  must  be  expanded  in  the  X and  Z 
directions  in  order  to  effect  containment,  but  the  depth  of 
four  layers  seems  sufficient  to  contain  all  sputtering 
events  at  energies  up  to  1 keV. 


D.  COMPARISON  OF  LATTICE  ORIENTATIONS 


comparison  of  ths  sputtering  ratio's  obtained  at  each 
orientation  vs.  ETHR  for  the  600-eV  and  1-keV  runs  yielded 
additional  insight  into  a possible  energy  range  for  the 
surface  binding  energy.  Examination  of  figures  3,5,  and  7 
shows  the  (100)  orientation  to  be  the  densest  orientation 
with  highest  probability  of  momentum  reversal.  The  (111) 
orientation  has  the  next  highest  probability,  and  the  (110) 
lattice  is  the  least  dense  with  high  probability  of  deep 
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penetration.  This  observation  leads  to  the  conclusion 
the  sputtering  ratios  should  decrease  in  the  order:  (100) , 

(111)/  (HO).  Figures  40  and  41  show  that  the  sputtering 

ratios  of  the  (111)  and  (110)  orientations  reversed  their 
order  of  sputtering  yields  at  a value  of  4 eV  in  both  the 
600-eV  and  1-keV  simulations.  The  600-eV  simulation  showed 
another  reversal  of  order  at  a binding  energy  of  9 e'7. 

These  data  points  would  then  indicate  that  the  most  probable 
range  of  binding  energy  exists  between  4 eV  and  9 eV. 

E.  MOMENTUM  TRANSFER  IN  THE  (100)  ORIENTATION 

Simulation  programs  were  run  with  bombarding  energies  of 
600  eV  and  1 keV  for  both  8x4x8  fee  copper  and  8x4x8  silicon 
diamond  lattices  with  the  impact  point  displaced  0.3  LU  in 
the  X and  Z directions  from  the  impact  atom's  center. 
Single  timestep  listings  of  the  atoms  were  obtained,  and  the 
spread  of  momentum  through  the  crystals  was  traced  by 
following  the  spread  of  lattice  atoms  which  possessed  total 
energies  above  0.1  eV.  The  traces  were  terminated  when  a 
boundary  atom  was  found  to  possess  a total  energy  greater 
than  0.1  eV.  In  addition,  simulations  were  run  in  which  the 
lower  energy  spread  of  momentum  (total  energy  > 0.0)  in  the 
top  layer  of  both  orientations  was  traced.  Figure  42  shows 
the  numbering  scheme  of  the  (100)  fee  crystal. 

1 • ££Q.Z§JL  Results 
a.  Copper 

Figures  43-43  show  the  spread  of  momentum 
through  the  fee  lattice  obtained  in  the  600  eV  simulation. 
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The  momentum  was  seen  to  travel  along  the  close-packed  (110) 
chain  in  the  surface  layer,  and  correspondingly  along  the 
adjacent  (110)  chains  in  the  second  layer.  by  timestep  27 
(1.345  x 10-14  sec)  the  upper  boundary  atom  was  reached  in 
the  surface  layer.  All  atoms  moving  with  an  energy  greater 
than  0.1  eV  were  observed  to  be  boundary  atoms  of  the  (110) 
channel  formed  by  the  major  diagonal  in  the  crystal. 


b.  Silicon 


Figures  49-57  show  the  momentum  spread  in  the 
(100)  diamond  lattice.  The  momentum  spread  here  was  also 
observed  to  travel  predominantly  along  the  boundary  atoms  of 
the  (110)  channel  formed  along  the  major  diagonal.  However, 
as  atoms  96,102,114,  and  120  began  to  move,  they  transfered 
momentum  outward  from  the  (110)  channel  into  neighboring 
channels.  The  momentum  spread  was  observed  to  be  very 
orderly,  and  the  mechanism  of  the  spread  was  observed  to  be 
the  obvious  fact  that  a moving  atom  affects  the  motion  of 
its  nearest  neighbors.  The  outward  spread  of  momentum  from 
the  first  layer  to  the  second  was  initially  slowed  in  the 
third  layer,  since  only  atoms  on  the  major  (110)  diagonal 
commenced  moving  first.  This  effect  was  again  found  to  have 
an  obvious  cause,  since  atoms  on  the  major  diagonal  were 
being  affected  by  the  motion  of  two  nearest  neighbors 
simultaneously.  An  extension  of  the  figure  at  timestep  72 
into  the  fourth  layer  showed  that  motion  on  the  three 
diagonals  of  the  third  layer  had  caused  motion  on  the  two 
diagonals  in  the  fourth  layer,  resembling  very  closely  the 
motion  of  the  second-layer  atoms.  It  is  assumed  that  if  a 
fifth  layer  had  been  added,  the  first  atoms  to  move  in  that 
layer  would  have  been  on  the  major  diagonal,  resembling  the 
motion  of  the  first  layer.  The  boundary  was  reached  in  the 
second  layer  at  timestep  72  (3.549  x 10-14  sec). 
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The  transfer  of  momentum  along  the  (110)  chain  in  the 
surface  layer  was  broken  by  the  sputtering  of  atom  57. 


2 • 1 0 0 0- e V Results 


a.  Copper 

Figures  58-63  show  the  momentum  spread  in  the 
fee  (100)  crystal  at  a bombarding  energy  of  1 keV.  Again, 
momentum  information  was  observed  to  travel  in  the  surface 
layer  (110)  chain;  however,  the  second  and  third-layer 
transfer  patterns  showed  an  irregularity  not  found  in  the 
600  eV  simulation.  This  irregularity  was  found  to  result 
from  the  passage  of  the  Ar  ion  between  atoms  53  and  57, 
forcing  these  atoms  in  opposite  directions.  This  collision 
started  a momentum  transfer  along  a (TlO)  chain  just  as  the 
initial  impact  started  a transfer  along  the  surface  (110) 
chain.  The  boundary  was  reached  at  timestep  29  (1.196  x 
IQ-14  sec) . 

b.  Silicon 

Figures  64-73  illustrate  the  momentum  transfer 
in  the  diamond  lattice  at  1 keV  bombarding  energy.  The 
mechansira  progressed  in  a manner  similar  to  the  600  eV 
simulation  with  the  exception  of  two  off-diagonal  atoms 
moving  in  the  surface  layer.  It  was  found  that  atom  42  was 
being  moved  by  the  slight  upward  movement  of  atom  113,  and 
atom  56  was  being  moved  by  the  slight  upward  movement  of 
atem  119.  This  then  indicated  that  the  momentum  is 
eventually  spread  in  all  directions  by  nearest-neighbor 
interactions  in  the  basic  tetranedra.  Once  again  the 
momentum  transfer  in  the  surface  (110)  chain  was  broken  by 
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the  sputtering  of  atom  57.  The  boundary  was  reached  at 
timestep  87  (3.337  x 10~14  sec). 

A comparison  of  the  momentum  transfer  in  the 
(100)  fee  crystal  and  the  (100)  diamond  lattice  showed  that 
the  spread  of  momentum  in  the  diamond  lattice  was  typically 
deeper,  and  that  the  information  was  spread  in  orthoganol 
directions  because  of  the  orientation  of  the  silicon  bonds 
in  the  basic  tetrahedra.  The  collision  chains  in  the 
diamond  lattice  seemed  to  be  more  easily  terminated  than  in 
the  fee  lattice,  but  the  information  spread  was  still 
orderly,  and  could  not  be  represented  by  a random  cascade  of 
collisions  such  as  would  be  seen  in  an  amorphous  solid. 

The  information  spread  to  the  lattice  boundaries 
was  2-3  times  slower  in  the  diamond  lattice  than  in  the 
corresponding  size  fee  lattice,  indicating  perhaps  that 
smaller  lattice  sizes  in  the  diamond  lattice  than  in  the  fee 
lattice  may  prove  to  sufficiently  contain  a sputtering 
event . 

3.  Precursor  Motion 

In  previous  simulation  work  with  fee  crystals,  it 
was  found  that  the  simulation  conditions  rapidly  transfer 
momentum  to  the  boundaries,  causing  very  slight 
displacements  of  all  atoms  in  the  crystal.  These  very  small 
displacements,  however,  affect  the  subsequent  motions  of  all 
atoms  because  of  the  program's  sensitivity  to  very  small 
displacements  and  velocities.  These  small  motions  would  be 
masked  by  the  thermal  motions  of  the  atoms  in  a real 
lattice,  but  they  are  one  of  the  irritating  problems  which 
must  be  addressed  in  the  simulation  programs.  In  order  to 
compare  the  effects  of  these  displacements  in  the  fee 
crystal  and  the  diamond  lattice,  a listing  of  all  atoms  in 
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the  top  layer  was  obtained  for  both  lattices,  and  the 
"precursor  motion"  of  all  surface  atoms  were  traced  until 
the  boundaries  of  the  surface  layers  were  reached. 

Figures  74  and  75  show  the  spread  of  momentum 
through  the  first  layer  of  an  8x4x8  fee  crystal  impacted  by 
a 600  eV  Ar  ion  at  a point  displaced  0.3  LU  in  the  X 
direction  and  0.3  LU  in  the  Z direction  from  the  center  of 
atom  21.  As  expected,  the  momentum  traveled  along  the  (110) 
close-packed  chain,  reaching  both  boundaries  within  five 
timesteps  (1.684  x 10-15  sec). 

Figures  76  and  77  show  the  advance  of  the  precursor 
motion  in  the  (100)  diamond  lattice  with  600  eV  bombarding 
energy.  It  was  observed  that  one  timestep  (approximately  4 
x 10  ~ie>  sec)  was  sufficient  time  to  transfer  momentum  to 
the  next-nearest  neighbor.  The  orthoganol  placement  of 
next-nearest  neighbors  in  the  tetrahedra  than  allowed  energy 
transfer  in  all  directions,  and  the  outward  spreading 
"ripple"  cf  momentum  reached  the  surface  boundaries  in  five 
timesteps  (1.894  x 10~15  sec).  The  diamond  lattice  and  the 
fee  lattice  exhibited  very  different  behaviors  in  these 
small-scale  energy  transfers,  but  the  speed  of  momentum 
transfer  to  the  boundaries  was  approximately  the  same. 
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IV. 


DISCUSSION 


A.  (100)  ORIENTATION 


The  (100)  orientation,  being  the  closest-packed  of  tne 
three  lattice  orientations  required  the  longest  computer  run 
times  and  possessed  the  highest  sputtering  ratios.  At  a 
bombarding  energy  of  600  eV,  the  sputtering  events  were 
confined  to  the  first  three  layers  (Fig.  16);  however,  nine 
atoms  located  on  t he  X or  Z boundaries  were  sputtered  in  the 
first  layer.  This  indicates  that  at  energies  as  low  as  600 
eV  an  8x4x8  crystal  size  is  inadequate.  The  1 keV 
simulation  sputtered  13  atoms  located  on  the  X or  Z 
boundaries  (Fig.  19),  11  in  the  first  layer,  two  in  the 
second  layer,  and  one  in  the  third  layer.  In  addition,  a 
low-energy  sputtering  evant  originated  in  the  fourth  layer. 
This  indicates  that  the  lattice  size  must  be  extended  in  the 
Y direction  if  higher  energy  bombardments  are  to  be 
performed. 

The  sputtering  mechanisms  observed  in  the  (100) 
orientation  (Fig.  18)  were  essentially  the  same  as  those 
observed  in  the  fee  copper  crystal.  The  one  significant 
exception  was  the  sputtering  of  atom  32  which  was  caused  by 
the  movement  of  a third-layer  atom  in  a direction 
essentially  parallel  to  the  X,Z  plane.  The  tetrahedral 
arrangement  of  the  silicon  atoms  allowed  momentum  transfer 
not  only  to  the  nearest  neighbor  (atom  32)  , but  also  to  the 
next-nearest  neighbors  (atoms  31  and  41)  . This  type  of 
interaction  demonstrates  the  fact  that  momentum  information 
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can  be  efficiently  transferred  from  the  surface  into  deeper 
crystal  layers,  and  conversely  from  deeper  crystal  layers  to 
the  surface,  allowing  these  deep  layers  to  have  a major 
effect  on  surface  layer  interactions.  This  type  of 
interaction  is  another  indication  that  the  crystal  depth  is 
inadequate. 

Examination  of  the  sputtering  ratios  obtained  at  600  eV 
and  1 keV  (Figs.  17  and  20)  showed  a rapid  increase  in 
low-energy  sputtering  events  below  a surface  binding  energy 
of  3 eY.  The  fourth-layer  sputtering  event  observed  at  1 
keV  bombarding  energy  was  below  the  3 eV  threshold,  and 
therefore  is  not  considered  a significant  sputtering  event. 

B.  (110)  ORIENTATION 

The  (110)  orientation  exhibited  the  same  containment 
problems  as  the  (100)  orientation.  At  600  eV,  sputtering 
events  were  confined  to  the  first  two  layers  (Fig.  23)  , but 
four  atoms  on  the  X or  Z boundaries  were  sputtered.  At  1 
keV,  five  atoms  were  sputtered  on  the  X or  Z boundaries 
(Fig.  26),  but  all  sputtering  events  still  occured  in  the 
first  two  layers.  With  the  bombarding  energy  set  at  2 keV, 
only  two  X and  Z boundary  atoms  were  sputtered  (Fig.  29), 
but  a fourth-layer  sputtering  event  occured  with  an  ejection 
energy  of  2.88  eV.  This  indicates  that  the  number  of  layers 
used  in  the  (110)  simulations  must  be  increased  at 
bombarding  energies  above  1 keV.  In  addition,  the  lattice 
dimensions  must  be  increased  in  the  X and  Z directions, 
especially  at  lower  energies. 

The  fact  that  the  momentum  in  the  (110)  orientation  is 
directed  deeper  at  higher  energies  was  observed  in  the 
mechanism  traces.  The  mechanisms  were  again  the  same  as 
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those  observed  by  Harrison  et  al  [12],  with  the  exception  of 
the  skip  mechanism.  The  skip  mechanism  rapidly  transferred 
energy  rather  large  distances  across  the  lattice  surface  and 
indirectly  caused  additional  sputtering  events  by  imparting 
high  velocities  to  at  least  two  other  surface  atoms  (Figs. 
25  and  28)  . Without  the  relatively  open  spaces  created  by 
the  (110)  channels,  and  without  the  placement  of  the  atoms 
in  the  (TlO)  chains  located  in  the  lower  half-layer  of  the 
surface,  this  mechanism  would  not  be  possible.  The  surface 
demonstrated  a rather  "brittle”  quality,  in  that  at  lower' 
energies  the  "skipping"  atoms  were  reflected,  but  at  2-keV 
bombarding  energy,  the  "skipping"  atom  was  deflected  deeply 
into  the  crystal.  This  deep  deflection  eliminated  the 
surface  sputtering  of  the  "skipping"  atom,  and  caused  high 
energy  interactions  in  the  lower  layers  (Fig.  31)  . 


The  comparison  of  the  sputtering  ratios  obtained  at  600 
eV,  1 keV,  and  2 keV  plotted  vs.  ETHR  (Fig.  32)  showed  a 
cross-over  of  the  sputtering  ratio  curves  at  binding 
energies  of  3 and  7 eV . A possible  explanation  for  this 
type  of  behavior  may  exist,  however.  As  the  energy  of  the 
bombarding  ion  is  increased,  the  target  cross  sections  of 
the  lattice  atoms  increase  to  a point,  and  then  become 
smaller.  This  type  of  behavior  could  explain  the  increase 
in  sputtering  ratio  observed  when  the  energy  was  increased 
from  600  eV  to  1 keV,  and  could  also  explain  the  lower 
sputtering  ratio  obtained  at  2 keV.  In  Fig.  32,  the  region 
between  3 and  7 eV  shows  this  type  of  behavior. 

C.  (Ill)  ORIENTATION 


At  both  600  eV  and  1 keV  bombarding  energies,  the 
sputtering  events  in  the  (111)  orientation  were  confined  to 
the  first  three  layers  (Figs.  34  and  37) ; however,  nine  X or 
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Z boundary  atoms  were  sputtered  at  600  eV,  and  11  were 
sputtered  at  1 keV . The  location  of  the  atoms  at  the  bottom 
of  the  "pseudochannels"  aided  in  containment  at  these  low 
energies,  but  it  is  expected  that  at  higher  energies  tne 
increased  crystal  damage  and  decreased  target  cross  sections 
of  these  blocking  atoms  may  necessitate  crystal  expansion  in 
the  Y direction. 

The  mechanisms  observed  in  the  (111)  orientation  were 
all  of  the  type  described  in  Harrison  et  al  [12],  but  the 
one  instance  of  reflection  of  an  atom  in  the  mechanism  trace 
(Fig.  36)  indicates  that  the  reflection  mechanism  is  more 
common  in  the  diamond  lattice  than  in  a fee  crystal.  It  is 
expected  that  at  higher  energies  this  mechanism  will  become 
less  common  because  of  lower-layer  damage. 

The  plots  of  sputtering  ratio  vs.  ETHR  in  the  (111) 
lattice  simulations  (Figs.  35  and  38)  show  the  rapid  rise  of 
low-energy  sputtering  events  below  3 eV  also.  This  type  of 
behavior  was  common  to  all  three  orientations,  and  seems  to 
be  a very  good  indication  that  the  binding  energy  is  above 
this  value. 


D.  SURFACE  BINDING  ENERGY 

Figures  40  and  41  show  the  comparison  of  the  sputtering 
ratios  of  the  three  lattice  orientations  plotted  vs.  ETHR. 
It  was  observed  that  above  a surface  binding  energy  of  4 eV, 
the  order  of  sputtering  ratios  at  both  bombarding  energies 
was  (100),  (111),  (110).  Considering  the  densities  of  the 
three  lattice  orientations,  this  was  the  order  expected  for 
the  sputtering  ratios.  In  the  600  eV  program,  tne 
sputtering  ratios  of  the  (110)  and  (111)  surfaces  were 
nearly  egual  above  a binding  energy  of  4 eV,  and  these  two 
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ratios  were  not  radically  different  than  the  sputtering 
ratio  of  the  (100)  orientation.  In  1963,  Southern,  Willis, 
and  Robinson  [23]  reported  their  investigation  of  the 
bombardment  of  Ge  and  Si  by  1-  to  5-  keV  Ar  ions.  Southern 
et  al  reported  that  they  did  not  find  any  significant 
differences  in  the  yields  of  the  (100),  (110),  and  (111)  Ge 
lattices  except  possibly  at  higher  energies  where  the 
sputtering  ratios  seemed  to  decrease  in  the  order  (100), 
(111),  (110),  with  the  (111)  and  (110)  yields  being  quite 
similar.  Figures  40  and  41  show  this  type  of  behavior  also. 
Considering  experimental  uncertainties  and  the  possible 
differences  in  the  sputtering  yields  of  Ge  and  Si,  it  can  be 
seen  that  the  sputtering  ratios  might  have  been  equal  within 
experimental  accuracy  at  lower  energies,  but  as  the 
bombarding  energy  was  increased  the  ordering  of  the 
different  orientations  became  more  evident.  The  interesting 
fact  is  that  by  adjusting  ETHR  to  a value  greater  than  4 eV, 
the  proper  ordering  of  the  sputtering  yields  is  obtained. 
Southern  et  al  also  reported  the  sputtering  ratio  of  the 
(111)  Si  crystal  orientation  at  1-keV  bombarding  energy  was 
approximately  0.45.  Figure  33  shows  that  this  sputtering 
ratio  was  obtained  in  the  1-keV  simulation  when  the  surface 
binding  energy  was  adjusted  to  a value  of  6-7  eV. 
Additionally,  Wehner  et  al  [2]  reported  the  sputtering  ratio 
of  polycrystalline  silicon  (uncorrected  for  secondary 
electron  emission)  to  be  approximately  0.52  at  a bombarding 
energy  of  600  eV.  The  mean  sputtering  ratios  of  the  three 
orientations  was  computed  for  the  600  eV  Ar  bombardment,  and 
was  found  to  be  0.567  at  4 eV  binding  energy  and  0.467  at  5 
eV.  These  rough  calculations  show  that  the  sputtering 
ratios  obtained  in  the  simulations  agree  reasonably  well 
with  the  available  experimental  data.  Anderson  [24] 
reported  obtaining  a sputtering  ratio  of  roughly  0.2  for 
the  Si  (110)  surface  bombarded  by  Ar  at  200  eV.  Although 
this  result  can  lead  to  no  definite  statement  about  the 
simulation,  it  is  reassuring  that  this  sputtering  ratio  is 
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below  the  ratios  obtained  in  the  (110)  surface  600  eV 
simulation  for  binding  energies  between  1 eV  and  10  eV.  A 
reasonable  assumption  then,  based  on  the  proper  ordering  of 
the  (100),  (111),  and  (110)  sputtering  ratios  between  4 eV 
and  9 eV,  and  the  reasonable  ordering  of  the  (110)  ratios 
obtained  at  600  eV,  1 keV,  and  2 keV  between  3 eV  and  7 eV, 
is  that  the  surface  binding  energy  is  between  4 eV  and  7 eV. 


E.  MOMENTUM  TRANSFER  STUDIES 

The  study  of  the  spread  of  energies  greater  than  0. 1 eV 
in  the  (100)  diamond  lattice  showed  that  the  momentum  was 
initially  transferred  along  the  (110)  chains,  but  spread 
outward  from  these  chains  because  of  the  placement  of 
nearest  neighbors  on  the  tetrahedron  diagonals.  The  lattice 
boundaries  were  first  reached  in  the  lower  layers  in  the 
diamond  lattice  (Figs.  57  and  73)  , as  compared  to  the 
surface  layer  in  fee  copper  (Figs.  48  and  63).  The  momentum 
spread  to  the  boundaries  reguired  2-3  times  more  timesteps 
in  Si  than  it  did  in  Cu.  This  indicates  then  that  the 
silicon  lattice  may  contain  all  sputtering  events  utilizing 
a smaller  lattice  than  the  fee  case  at  corresponding 
bombardment  energies.  Since  the  boundaries  are  reached 
later,  these  boundary  conditions  will  have  a smaller  effect 
upon  the  sputtered  atoms,  and  since  the  boundaries  are 
reached  in  lower  layers,  they  will  have  less  effect  on  the 
predominant  surface  sputtering  events. 

The  "precursor  motion"  in  the  silicon  lattice  (Figs.  76 
and  77)  was  observed  to  spread  rapidly  in  all  directions, 
rather  than  along  a (110)  chain  as  it  did  in  copper  (Figs. 
74  and  75) . It  was  originally  expected  that  momentum  would 
be  transmitted  deeper  into  the  lattice,  but  would  spread 
outward  to  the  X and  Z boundaries  more  slowly.  The  fact 
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that  every  atom  in  the  surface  layer  was  displaced  by  the 
fifth  timestep  was  surprising  but,  in  retrospect,  should 
have  been  expected.  These  small  displacements  and  atom 
vibrations  are  masked  by  thermal  effects,  since  at  the  end 
of  the  fifth  timestep  the  largest  kinetic  energy  possessed 
by  any  atom  other  than  the  impact  atom  was  on  the  order  of 
10-3  eV,  and  atomic  displacements  were  on  the  order  of  10~* 
LU  (2.72  x 10-4  angstroms) . This  computer  artifact  then  has 
no  physical  significance,  but  it  does  affect  the  subsequent 
motions  of  all  atoms  in  the  lattice.  These  small 
displacements  cause  the  atoms  to  oscillate  about  their 
equilibrium  positions  until  a higher-energy  interaction 
causes  a larger  displacement.  The  period  of  these 
oscillations  ( 10-16  sec)  is  below  that  which  would  be 
expected  in  thermal  vibrations.  This  motion  is  merely 
numerical  noise,  and  must  be  tolerated,  but  not  ignored  when 
constructing  a simulation  model.  The  observed  precursor 
motion  does,  however,  demonstrate  that  different  lattice 
configurations  do  transfer  momentum  information  in  radically 
different  ways. 


V. 


CONCLUSIONS  AND  P2C0MM BN DATIONS 


The  preliminary  development  and  investigation  of  the 
silicon  lattice  produced  physically  reasonable  results 
utilizing  an  acceptable  amount  of  computer  time.  The 
simulations  performed  with  ion  bombarding  energies  of  600 
eV,  1 keV,  and  2 keV  resulted  in  the  identification  of  the 
skip  sputtering  mechanism,  which  was  not  identified  by 
Harrison  et  al  [12]  in  the  fee  copper  lattice.  The 
reflection  mechanism  was  found  to  be  more  common  in  the 
diamond  (111)  orientation  than  in  the  closely-packed  fee 
crystal.  Adjustment  of  the  binding  energy  to  a value 
between  4 and  7 eV  resulted  in  the  proper  ordering  of  the 
sputtering  ratios  of  the  (100),  (110),  and  (111)  lattice 
orientations.  The  sputtering  ratio  of  the  (110)  lattice  was 
found  to  decrease  as  the  bombarding  energy  was  increased  to 
2 keV  as  a result  of  the  decreased  target  cross  sections  of 
the  lattice  atom.  A lattice  size  of  256  atoms  was  found 
inadequate  to  contain  all  sputtering  events  at  a bombarding 
energy  as  low  as  600  eV,  but  was  deemed  large  enough  to 
conduct  preliminary  bombardment  simulations.  The  spread  of 
momentum  in  the  diamond  lattice  was  found  to  be  2-3  times 
slower  in  the  diamond  lattice  than  in  the  fee  lattice.  The 
initial  studies  performed  should  be  continued  and  expanded 
upon.  Wehner  et  al  [1]  have  gathered  data  on  the 
temperature  dependence  of  the  sputtering  yields  of  silicon, 
and  the  ejection  patterns  obtained  from  the  (100),  (110), 
and  (111)  orientations.  The  computer  program  presently  has 
the  capability  to  randomly  displace  the  lattice  atoms  as  a 
function  of  lattice  temperature,  and  with  slight 
modifications  could  be  made  to  gather  statistical  data  on 
the  ejection  directions  of  the  sputtered  atoms.  Since 
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experimental  data  are  readily  available  in  these  areas,  the 
extension  of  the  silicon  investigation  into  these  areas  is 
recommended . 


A.  LATTICE  MODEL 


The  lattice  models  of  the  three  orientations  do  not 
presently  include  the  bulk  relaxation  of  the  lattice  surface 
layers.  Harrison  et  al  [14]  included  surface  layer 
relaxation  in  their  improved  fee  copper  model  and  found  that 
emission  directions  were  modified  slightly  by  the  relaxation 
process,  and  the  energy  distribution  of  the  sputtered  atoms 
was  modified  slightly.  it  is  recommended  that  surface 
relaxation  be  applied  to  the  silicon  lattice  also. 

The  potential  approximations  utilized  in  the  simulation 
programs  yielded  physically  reasonable  lattice  behaviors  at 
the  bombardment  energies  investigated;  however,  the 
relative  error  of  the  Si-Si*1  potential  increases  above  ten 
percent  at  potential  energies  above  10  keV.  If  high  energy 
bombardments  are  to  be  simulated,  it  would  be  advisable  to 
improve  the  approximation  at  small  atomic  separations. 
Inclusion  of  an  attractive  potential  region  could  prove  to 
be  extremely  difficult,  and  perhaps  over sophisticated.  In 
order  to  properly  include  an  attractive  potential  region, 
the  directional  effect  of  the  potential  function  would  also 
have  to  be  included.  If,  in  fact,  this  could  be 
accomplished  it  is  expected  that  the  computer  run  times 
involved  in  simulating  the  bombardment  of  a suitable  lattice 
size  could  become  prohibitive.  It  is  expected  that  the 
inclusion  of  an  attractive  potential  region  would  not 
materially  affect  the  identification  of  the  sputtering 
mechanisms,  but,  as  pointed  out  by  Harrison  et  al  [14],  the 
inclusion  of  an  attractive  potential  would  eliminate  the  use 
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of  a binding  energy  term,  which  is  a convenient 
parameter. 

The  problem  of  adequate  containment  of  all  sputtering 
events  remains  unresolved.  While  the  256  atom  lattice  was 
deemed  large  enough  for  preliminary  studies,  it  is 
recomended  that  the  lattice  size  be  increased  until  farther 
lattice  size  increases  do  not  affect  the  sputtering  ratio 
obtained  at  a given  bombarding  energy.  The  value  of  ETHH 
can  then  be  adjusted  to  obtain  a sputtering  ratio  which 
conforms  to  available  experimental  data. 


Preliminary  diamond  lattice  investigations  snowed  that 
reasonable  sputtering  ratios  wer°  obtained,  and  correct 
ordering  of  the  sputtering  ratios  of  the  (100),  (110),  and 
(111)  orientations  were  obtained  with  the  binding  energy  set 
between  4 and  7 eV . Further  investigations  of  sputtering 
with  larger  lattice  sizes  should  be  performed,  and  the 
results  correlated  with  available  experimental  data  in  order 
to  further  localize  the  binding  energy  value. 
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The  (111)  orientation  showed  the  reflection  of  a surface 
impact  atom  from  second  and  third- layer  atoms,  a mechanism 
reported  only  rarely  in  the  fee  lattice.  In  all 
orientations,  the  efficient  transfer  of  momentum  to  all 
atoms  in  the  basic  silicon  tetrahedron  was  noted.  It  is 
recommended  that  the  energy  dependence  of  the  skip  and 
reflection  mechanisms  be  investigated  in  order  to  verify 
their  decreasing  importance  at  higher  bombarding  energies. 

The  sputtering  ratios  of  all  three  orientations,  but 
especially  the  (110)  orientation,  should  be  studied  as  a 
function  of  bombarding  energy.  The  point  at  which  the 
decreased  target  cross  sections  of  the  lattice  atoms  cause 
decreased  sputtering  ratios  should  be  determined,  and  these 
data  correlated  with  available  experimental  data.  This  will 
aid  in  the  determination  of  the  surface  binding  energy. 

The  investigation  of  the  momentum  spreading  mechanism  in 
the  (100)  orientation  showed  that  the  momentum  initially 
spreads  along  the  (110)  chains,  much  as  it  does  in  the  fee 
lattice,  but  the  tetrahedral  arrangement  of  the  lattice 
atoms  spreads  the  momentum  orthogonal  to  the  (110)  diagonal, 
and  allows  the  efficient  transfer  of  momentum  to  the  lower 
layers.  3oundary  atoms  were  reached  in  lower  layers  before 
the  boundary  atoms  in  the  surface  layer  were  reached.  The 
momentum  spread  required  2-3  times  as  many  timesteps  in  the 
8x4x8  diamond  lattice  than  it  did  in  an  8x4x8  fee  copper 
lattice  impacted  at  the  same  energy.  The  small  amplitude 
"precursor  motions"  in  the  surface  layer  were  found  to 
spread  to  all  directions  in  the  surface  layer,  rather  than 
along  the  (110)  diagonal  as  it  did  in  the  fee  lattice.  The 
speed  at  which  the  "precursor  motions"  spread  to  the 
boundaries  was  found  to  be  approximately  the  same  in  the 
diamond  lattice  and  the  fee  lattice. 
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Further  investigation  of  the  momentum  spread  in  the  (110) 
and  (111)  may  provide  additional  insight  into  the  dependence 
of  the  sputtering  ratios  upon  the  atomic  arrangement  of 
various  lattice  types. 

Further  investigation  of  the  diamond  lattice  could  prove 
fruitful  in  the  understanding  of  the  problems  and  mechanisms 
involved  in  semiconductor  etching,  ion  implantation,  and 
depth  profiling. 
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APPENDIX  A 


DIAMOND  LATTICE  SUBROUTINES 


COT 

Direction  cosine  of  the  bombarding  ion's  velocity. 

DBIHAX,  DBZMAX 

Upper  boundary  limits  for  the  target  impact  area. 

DX,DX , DZ 

Unsealed  displacement  distance  components  of 
second  atom  at  a given  lattice  site. 

I,J,K 

Integer  index  variables. 

ILL 

The  number  of  the  last  atom  in  the  first  layer. 

IN 

Index  integer  variable  in  DL111  subroutine  atom 
placement. 


IT#JT,KT 

Index  variables  used  in  the  generation  of  an 
atom's  X , Y , Z coordinates. 

ITT 

An  odd-even  integer  used  to  determine  atom  site 
location . 
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IX,  IT,  IZ 

The  number  of  x,y,z,  planes  in  the  crystal. 


JTS 

Variable  use!  in  DL111  subroutine  to  determine  the 
Z coordinate  of  a given  atom. 

JTT 

Index  variable  in  DL111  subroutine. 

L 

Number  assigned  to  the  second  atom  at  a given 
lattice  site. 

LD 

Number  of  the  next-to-highest  numbered  mobile 
atom. 

LL 

The  highest  numbered  atom  in  the  entire  crystal. 

H 

Number  assigned  to  the  first  atom  at  a given 
lattice  site. 

HI 

Distance  from  lattice  impact  point  to  the 

bombarding  ion's  initial  position. 

BOE 

Lattice  nearest-neighbor  distance. 

BXI  (I)  , BTI  (I)  , RZI  (I) 

Initial  position  of  the  Ith  atom  in  the  lattice 
(used  to  reposition  the  atom  before  the  next 
impact)  . 
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used  to  convert 


SCI,  SCI,  scz 

The  x,y,z  lattice  scale  factors 
the  generated  lattice  to  match  a real  crystal. 


SI 

Reciprocal  of  the  square  root  of  2.0  (used  to 
avoid  repeated  division)  . 

SPDI,  SPDZ 

The  fraction  of  the  target  impact  area  that  the 
impact  points  are  incremented. 

SSPDZ 

SPDZ  divided  by  3.0  (used  in  DL 1 1 1 subro utine  to 
avoid  repeated  division)  . 

sscz 

SCZ  divided  by  3.0  (used  in  DL111  subroutine  to 
avoid  repeated  division)  . 

ZP 

Unsealed  z displacement  variable  used  in  DL111 
subroutine  to  center  the  lower-layer  atoms  in  the 
equilateral  triangles  formed  by  the  surface  atoms. 
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0L100  G c N E R A T E S A 01  A MONO  LATtICt  im  tH-  (iqO) 
ORIENTATION.  A SIMPLE  ^ AC  E CENTERED  C'JBIC  L'-TTTCE 
IS  GEN  RATED,  AND  THE  N AN  ADDITIONAL  A~OM  T ' ADDED  TO 
EACH  LATTICE  SITE  DISPLACED  ( 1 / 2 , I /2 , 1/ 2 , ) LU. 


COMMON/ C OMl/RXl 600) , P Y ( 600 ) *5  Z (600) * LC;JT (600) • LL , LD 
COMMON/ CnM2 /R 0 -, ROc2, ROE M, AC  , PAC , P?TC , PTC , PFRT' , c dtC , 

1EM,PFIV ,TRQT 

COMMON/ CDMA /IX, I Y, I Z , I XP , I Y P , I Z P , 5CX , SC Y , SC Z , I CE E D , 
2°EX,PBZ 

CCMMCN/C0M7/R1 ,LSS,SRX,SPZ,CnX,COY»COZ, YlAX (20), ILL 
COMMON  / CO M9  / IOX,  I OY,  ItZ,OXT,~»YT,DZ~,TPOX,tPCY,_POZ 
CCMM0N/C0M1  I/RX I ( 600)  ,RYI1600),RZI  (600) 

C CM MON/ C 0M2  3/  S PD X , SP  D Z , D 3X  M AX , D B ZMA  X 

SET  I V P A C T’  AR-A  PARAMETERS 

Lp  = 0 

SDCX=0.  125 
SPCZ  = 0.  125 
CRXMAX=  1.0 
C E Z MA  X = 0 . 0 
SCX=1.0 
SCY  = 1 .0 
SCZ  = 1 .0 


BEGIN  LATttc=  GENERATION 

THC  L AT  T T C E IS  GENERATED  FIRST  IN  THE  X DIP7CTI0N, 

T FEN  IN  ~ HE  Z DIRECTION.  U00:  COM  PL  _T0N  °F  A L*Y-R, 
A NCW  LAYS'?  IS  BEGUN  IN  TH - INCREASING  Y 0IpcC~7QN. 


M = 2 
M = 3 

IF(LSS.EQ.3)  M = 3 
DX=0. 5 
D Y = 0 . 5 
DZ=0. 5 
JT  = 0 

Y =-SC Y 

DC  60  J =1  , I Y 

Y = Y+SCY 
K T =0 

z=-srz 

DC  50  K =1 , T Z 
Z=Z+SCZ 
IT  = 0 
X=-SCX 

00  AO  1=1, TX 

x=x+scx 

ITT  = I t+Jt  + kt 

I F ( ITT- ( ITT/2) *2 ) 39,30,39 
30  PX(M)=X 
RY(M) =Y 
R Z ( M ) =Z 
RXI ( M ) =X 
RYI (M )=Y 
P Z I ( M ) = Z 
L=M 
M =M  + 1 

P X ( M ) =P  X ( L ) +0X 
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PY  (M)  =R  Y{  L ) +DY 
RZ(M)=PZ(l)+DZ 
P XI ( M ) =RX  C M ) 

FYI  ( M ) = RY  ( M ) 

PZI  ( M ) =RZ  ( M ) 

M =M  + 1 

39  C1DNT  I NU  E 
I T = I T+  1 

40  CONTINUE 
KT=KT+1 

5C  CONTINUE 

IF(JT.-g.o)  I LL  = 1 

J~=JT  + 1 

Ic< JT.=Q.IDEEP)  LD=m-1 
60  CONTINUE 
LL=M-I 

Ic(LD.“G.O  ) LD=LL-1 
IF(CCY.EO.O.O)  GO  -D  120 
P 1=RC  E/AB  5 ( COY) 

RETURN 
120  R 1 =1 . 5 
RETURN 
END 


SLeROUTINc  CL  U 0 


CLUO  G E’'!  Ec  A“cc  a DIAMOND  L A T T I C E IN  TH-  (110) 
ORIENTATION.  A CHAIN  OF  ATOMS  IS  Ftqjt  GENERA"7  0 IN 
T HE  X DIRECTION,  AND  tHcN  tH-  SAME  CHA7N  IS 
GENERATED  A FT  ER  BEING  DISPLACED  IN  TH:  Z DIRcTTTON. 
WHEN  A LAYER  IS  c I L L E D , TH-7  SAME  CHAINS  ARE  GENfRAT 
IN  THE  NEXT  LAY"R  STARTING  AT  DIFFERENT  I N I T I A L 
X AND  Z POSItthms 


C C MpMCN/  COM!  / R X ( 6 0 0 ) , RY(600  ) , = Z ( 600 ) , LCUT ( 60 0 ) , LL , L 0 
CCMM0N/C0M2/R0E  ,c0-2,R0EM,AC,:>AC  , PPTC  , P~  C , P P°~r  , -PT 

1FN,PFIV,TPQT 

COMMON/ COMA/ IX, IY, IZ, IX3, IYP, I Z3,SCX,SCY,SCZ,  I DEEP, 
2 F 8 X » ?,  B Z 

CCMMON/OOM7/R1,LSS,SPX,3PZ,COX,i"oy,COZ,  YLAX(20)  , I LL 
CCMMON/COM9/  IDX,  IDY,  I DZ , DXT  , DY  t , QZt  , T P 0 X , T POY , "P  "1 Z 
COMMON /COM l 1/PXT  (600)  ,RYI(60)),RZI(600) 

CCMMCN/ COM2  3/  S PDX , SPOZ, DBXM AX , D8ZMA X 
SET  IMPACT  AREA  PARAMETERS 

S I =1 . 0/ SQRT  (2.0) 

L C = 0 

SCX  = S I 

SC Y=S  I 

SCZ=1.0 

r3XMAX=1.0 

C 8 ZMA  X = 2 . 0 

SPDX=OPXMAX/5.0 

SP0Z  = DBZMAX/5 .0 

D X = S I 

CY  = 0.  3 

DZ=0.5 


E-C-IN  LATTICE  gcneration 


M=  3 

JT=0 

Y=-SCY 

DC  60  J = 1 , 1 Y 
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Y =Y+SCY 
K T =0 
2 = -SC  Z 

DO  50  K = 1 , I z 

z=z+scz 

17=0 

x=-scx 


DC  90  I =1 , I X 
X=X+SCX 

IR( I T— ( !T/2)*2) 

21,11,21 

11 

IF ( JT-( JT/2  )*2) 

90, 12,90 

12 

IF (KT-( KT/2)-2) 

90,30,90 

21 

IF(JT-( jT/2)*2) 

22,90,22 

22 

IF  (KT-(KT/2)*2) 

30,90,30 

3C 

R X ( M ) =X 
R X I ( M ) = X 
R Y ( M ) =Y 
RYI ( M ) = Y 
R Z ( M ) =Z 
F Z I ( M ) = Z 
L =M 

M=M  + l 

0 X ( M ) =R  X ( L ) +-0X 

R X I ( M ) = R X ( M ) 

RY  ( M)  =RY(  L ) +DY 

R Y I (M)=RY( M ) 
RZ(M)«?Z(L)+nZ 
PZI  ( M ) =P  Z { M ) 
M=M  + 1 

90 

17  = 17+1 
K T =KT  + 1 

53 

CONTINUE 
IF(JT.=g.o)  ILL 
J7=JT+1 

■=  M - 1 

I F( JT.RQ.IDERD) 

L0=M-1 

60 

CONTINUE 

L L =V-  1 

IF(LD.EC.O)  ld= 

LL-l 

I F ( COY. £0.0.0) 

R 1=RQE/ A3S ( CnY) 
R CTUR  N 

GO  TO  120 

120 

R 1 = 1. 5 
P c TURN 
ENC 

SLPR0UTINE  DL111 


0L111  GENERATES  A DIAMOND  L A~7  ICE  TN  THE  (111) 

OR  I ENTA  T I ON  . THC  ATOMS  A R c PLACED  IM  PA  I 0 S ~ 7 "ACH 
LATTICE  SITE,  WITH  GEM=RA~IOM  p ! RST  IN  tHE  X AMD  Trl-M 
IN  THE  Z DIR£CtION.  WHEN  THE  cIRST  LAYER  tS  complet- 
THE  SECOND  LA  Yc  R ATOMS  ARE  PLACCD  OIREC'LV  H = N"ATH 
the  CENTROIDS  Oc  TH-  EOUILAT-ral  tRIANGLES  cr>RM  = n 
BY  THE  FIRST  LAY=R  ATOMS.  ~HC  THIR0  LAY^R  ATOMS  ARE 
PLACED  B ENE  ATH  THE  S -COND  L AYcR  CENTROIDS.  "HE  =0UR7H 
LAYER  IS  A REPEAT  OF  THE  FIRST  LAYER. 


CCMMON/CGMI  /RX(  60  0)  . R Y (600 ) , R Z ( 600 ) , LC'JT  ( 600  ) , LL  , LD 
CCMMDN/C0M2/P  Oc,  R0F2,R0=M,  AC,  PAC  , PP  ~C  , P TC  , P FP  Tr  , F PTC  , 
1FM, PFIV .TPQT 

COMMON/COM  A /IX,  I Y , I Z , I X? , T Y P , I Z R , S CX , SC  Y , SC  Z , IDE  c P , 
2FBX,RBZ 

CCMMCN/C0M7/R1,LSS,SPX,$PZ,C0X,CDY,C0Z, YLAXI20) , ILL 
CCMM0N/C0M9/  IDX , IDY ,IDZ,DXTt DY’.OZT, TDOX, TPOY, TPOZ 
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CCMMCN/C0M11/PXI  (600)  , RYI  (600)  * Q Z I (600) 
CCMMCWCOM2  3/  S^DX , SPOZ , D8XMAX  , 08  ZM-\  X 
CCMMCN/C0M26/SSPQZ 

SET  impact  area  parameters 


ID  =0 

SCX=I. J/SQRT( 2. 0 ) 
SCY=2.0/SQ°T(3.0 ) 

SCZ  = SQR  T ( 1.5) 

SSC  Z=SC  Z/ 3 . 0 
DEZMAX=2.0 
C E XMA  X=  1 . 0 
SPDX=0BXM AX/5.0 
SPDZ=0BZMAX/10.0 
SSPDZ=SPDZ/3. 0 
DX  =0 . 0 
D Y=Rf)E 
CZ  = 0.0 

BCGIN  LA^TIC=  GENERATION 

JT  =0 
M = 2 

Y =-SC Y 

DC  60  J =1 , I Y 
Y= Y+SCY 
KT=0 

J'rS=JT+JT/3 

Z=-SCZ 

DC  50  K = 1 » I Z 
Z=Z+SCZ 

I T = J 

X=-SCX 

CC  40  1=1,  IX 

X=X+SCX 

IN=IT+JTS+KT 

I F ( IN— ( IN/2  )*2 ) 39,30,39 

30  RX(M)=X 

RXI  ( M ) = R X ( M ) 

PY(M) =Y 
PYI  ( M ) = R Y ( M ) 

IF( JT-3*( JT/3) ) 31,35,31 

31  JTT  =JT 

22  JTT  =JT7-3 

IF(JTT)  33,35,32 
33  JTT -JT~ +3 
ZP=JTT 

RZ(M)=Z+Z°*SSCZ 
RZI  ( M ) = R Z ( M ) 

GG  to  37 
35  R Z ( M ) = Z 

RZI (M)=RZ(M) 

37  L=w 
M =M+1 

P X ( M ) =P  X ( L ) +0X 
PXI (M)  = RX( M ) 

RY(M)  =P  Y ( L ) +0Y 
RYI (M ) =RY( M ) 

R Z ( M ) =R  Z ( L ) +0Z 
RZ!(M)=RZ(M) 

M =M+  1 

39  CONTINUE 

I T = I T+ 1 

40  CONTINUE 
KT=KT+1 

50  CCNTINUE 

IF(JTT.EQ.O)  I L L = M-  l 
JT=JT+1 

IF(JT .cO.IOE-P)  L 0=M- 1 
60  CCNTINUE 
L L =M- l 


8 1 


M7) 


/ 


APPENDIX  3 


ENERGY  SUBROUTINE 


A 

Exponential  Bohr  potential  parameter  in  Si-Si 

interaction . 

B 

Exponential  Bohr  potential  parameter  in  Si-Si 

interaction. 

DIST 

Variable  used  alternately  as  the  square  of  the 
distance  between  two  atoms,  and  as  the  distance 
between  two  atoms. 

DRX,DRY,DRZ 

X,Y,Z  components  of  the  displacement  vector 
between  two  atoms. 


EXA, EXB 

Exponential  Born-Mayer  potential  parameters  in 
Si-Si  interaction. 


I 


Index. 


IP 


Index. 
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J 


Index. 

LCOT  (I) 

Integer  which  indicates  whether  the  Ith  atom  has 
exited  the  sides  or  bottom  of  the  crystal,  has 
risen  from  the  crystal  surface,  or  is  still  within 
the  crystal. 

LD 

Number  of  the  next-to-highest  numbered  mobile 
atom . 

LL 

Highest  numbered  atom  in  the  lattice. 

PA , PB 

Bohr  exponential  coefficients  in  Ar-Si 
interaction. 

PAD, PBD 

Transition  potential  exponential  coefficients  in 
Ar-Si  interaction. 


PEXA , PEX3 

Born-Mayer  exponential  coefficients  in  Ar-Si 
interaction . 


PFIV 

0.5 

POT 

Potential  energy  generated  between  atoms  by  their 
interaction . 
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PPE(I) 

Potential  energy  of  the  Ith  atom. 

PPTC 

Ar-Si  interaction  potential  evaluated  at 

nearest- ne ighbor  separation. 

PTC 

Si-Si  interaction  potential  evaluated  at 

nearest-neighbor  separation. 

PV3M 

Born-!layer  term  of  the  Ar-Si  interaction 

potential. 

ROE 

N ear est- ne ighbor  distance. 

R0E2 

Nearest-neighbor  distance  squared. 

ROEA 

Cross-over  point  of  the  Bohr  and  Born-layer 

portions  of  the  Si-Si  interaction  approximation. 

ROEB 

Cross-over  point  of  the  Bohr  and  Born-Mayer 

portions  of  the  Ar-Si  interaction  potential 

approximation. 

BOEC 

Truncation  distance  of  the  Ar-Si  interaction 
potential  (set  at  nearest-neighbor  distance)  . 
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B0EC2 

Ar-Si  truncation  distance  squared. 

SAVE 

One-half  of  the  potential  energy  generated  by  the 
interaction  of  two  atoms. 

HX(I)  ,RY(I)  ,RZ(I) 

X , Y,  Z coordinates  of  the  Ith  atom. 


TPOT 


Total  potential  energy  of  all  atoms. 


XHAX,YHAX,ZHAX 

Maximum  X,Y,Z  coordinates  of  any  atom  which  could 
interact  with  the  atom  being  evaluated. 


XMIN#YMINrZMIN 

Minimum  X,Y,Z  coordinates  of  any  atom  which  could 
interact  with  the  atom  being  evaluated. 
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SLBROU”  INE  ENEFGY 


THIS  SUBR  OUTT N:  CAl.CULAT=S  THE  MUTUAL  D0T5NT I SL 
ENERGIES  GF  INTERACTING  ATOMS.  THIS  SUBROUTINE  4nj 
ACAPTIPN  GF  THr  MUNDKUR  HI0H-SP5EC  V = RC  I IN  *IITH 
TWO-BRANCH  APPROXIMATIONS  OF  "rHE  P0T5NT  I A L S . 


COMMnN/COMl/RX(60O)  ,RY(600)  ,RZ  (600)  tLC'JT  (60  0)  ,LL,L0 
COMMON/ COM2 /R  Or » R 0 E2 » ROEM , A C , 0 AC  , PPTC  *DTC,PcPTr ,FPTr , 
1CP ,PFIV ,TP0T 

CCMMCN/C^M3/rXA ,CX3 ,FXA ,P5XA,PrX3,PFXA , THERM 
COMMON/ C0M5/F X ( 600 ),FY( 500)  ,FZ(600),PP:<  600) 
CCMM0N/C0M22/R05 A, ROCB,RO=C .Rn-C2.CP0,CP I , C &2 , CP 3 , C p0 , 
2CC1,CF2,CGD1,CG0  2,CGB1,CG32 
COMMON/ COM2 4/ PA, PB,sAD,P3D, A, 3, PFAD.CVED 

3CM9AR0ING  I0N-LAT~!C5  ATOM  : NT -R AC" I ON  COMMUTATION 


1=1 

I F ( LCUT (I ) ,N=. 3)  GO  T0  200 

XN  TM=  RX ( I ) -R05C 

XMAX=RX( I ) + RQEC 

YM IN  = RY ( I ) -ROEC 

YMAX  = RY ( I I + ROEC 

ZMIM=RZ ( I )-ROEO 

ZM Ax=  RZ ( T ) + R0EC 

OC  195  J =2  » LL 

I c (LCUT ( J ) . Nc  . 3 ) GO  TO  195 
Ip (RY(J ) .GE.YMAX  ) GO  T0  195 
!F(RY(J).LC.YMIN)  GO  TO  195 
IF(PX(J l.G-.XMAX)  GO  T0  195 
I F ( RX ( J ).L5.XMIM)  GO  To  195 
!F(RZ(J ).G5.ZMAX)  GO  TO  195 
I F ( RZ ( J J.LE.ZMIN)  GO  TO  195 
DR  X =RX ( J ) -RX ( I ) 

DRY=RY( J )— RY( I ) 

DRZ=RZ( J)-RZ( I ) 

D I ST  = DR  X**2+0RY^^2+D3  Z**2 
IF(OIST .GE.R0EC2 ) GO  TQ  195 
D!ST=SQRT (DIST) 

°veM  = EXO(P5XA  + ?5XBT-DTST)-pPTC 
I F ( D I ST .GT.R0F3)  GO  '0  175 
P0T=FXP (PA+P3*D ! ST) /OI ST+PVBM 
GO  T0  180 

175  PCT=EXP ( 0 A0+P8D*D 1 ST ) +PVBM 
180  TPr,T=TpQT  + P0T 
SAVE=pc I V * P QT 
PPEtl  ) = PPE  ( D+SAV5 
P°F(J)=PPE(J) hSAVE 
195  CCNT I NU c 
2C0  CONTINUE 

LATTICE  ATOM  INTERACTION  COMOJT A”  IQM. 
CC  303  1=2,  LD 

I f ( LC  UT ( I ) . NE  . 3 ) GO  ’0  300 
IP=I+1 

XMAX  = RX  ( I ) +-ROE 
XMI N=RX ( I ) -ROE 
Y MA  X=  RY ( I ) + ROE 
YM I N = RY ( I ) -ROE 
ZM A X=RZ ( I )+RHF 
ZM IN=  RZ ( T ) -D0  E 
DC  295  J=  I P , LL 
I F ( LCUT  ( J ) . NE  . 0 ) GO  T0  295 
IF (RY ( J ) .G=  .YMAX)  GO  T0  295 
IF  (RY<  J ) . LE.YMI  N)  GO  T~>  295 
IF (RX ( J ) .05  .XMAX ) GO  "0  295 
IF(ox(j).L5.XM!N)  GO  TO  295 
IF(FZ( J J.G5.ZMAX)  GO  to  295 
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L 


I 


I P (RZ ( J ) . LE .ZMI M)  GO  T0  295 
OR  X =R  X ( J ) -R  X ( ! ) 

OPY=RY( J)-RY(  I ) 

DRZ=RZ( J)-RZ( I ) 

D !ST=0RX**2+nRY**2+P3 Z**2 
I F ( DI ST . GE . ROE2 ) GO  tq  295 
DI ST=SQRT(D  1ST) 

POT=PXP (EXA+FXB*OIST)-PtC 
IF(OTST.GT.ROEA)  GO  ~1  289 
PCT=EXP (A+B*DI ST ) /OI ST+POT 
280  tpqt  = tpCt  + °Ct 
SAVE=PFI V*PGT 
PPE  ( I ) = PPE  ( D+SAVE 
P°E ( J ) =PPE ( J ) +SAVE 
295  CONTINUE 
3JJ  CONTINUE 

DC  310  1=1, LL 

I F ( LC UT ( I ) .EQ. 3)  GO  "0  31  3 
T POT  = TPQT  + Pr’-  ( I ) 

310  CONTINUE 

return 

END 
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APPENDIX  C 


STEP  SUBROUTINE 


A 

Exponential  Bohr  potential  parameter  in  Si-Si 

interact  ion. 

AC 

Si-Si  exponential  force  coefficient  obtained  by 
dif f erentation  of  the  Born-Sayer  potential  term. 

B 

Exponential  Bohr  potential  parameter  in  Si-Si 

interact  ion. 

CVED 

Ratio  used  to  convert  forces  to  KKS  units. 

DPP 

Difference  between  the  separation  distance  of  two 
atoms  and  force  truncation  distance  (ROEC-DIST) . 

DIST 

Variable  used  alternately  as  the  square  of  the 
distance  between  two  atoms,  and  as  the  distance 
between  two  atoms. 

DRX,DRY,DRZ 

X,Y,Z  components  of  the  displacement  vector 
between  two  atoms. 
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EXA,EXB 

Exponential  Born-Mayer  potential  parameters  in 
Si-Si  interaction. 

FA 

Directional  force  component  (force  multiplied  by 
direction  cosine) . 

FBH 

Force  term  obtained  by  dif ferentation  of  the 
Born-Mayer  term  in  the  Ar-Si  interaction. 

FH 

Minimum  force  given  to  an  atom,  atoms  having 
forces  exerted  on  them  which  are  less  than  FM  have 
their  forces  set  to  zero. 

FOD 

Force  divided  by  distance. 

FORCE 

Magnitude  of  the  repulsive  force  between  two 
atoms. 

FPTC 

Si-Si  force  evaluated  at  nearest-neighbor 
distance . 

FI  (I)  ,FT(I)  ,FZ(I) 

X, Y,Z  force  components  exerted  on  the  Ith  atom. 

FXA 

Si-Si  exponential  force  coefficient  obtained  by 
dif ferentation  of  the  Bohr  potential  term. 
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r 


i 


Index. 


IP 


J 


Index. 


Index. 


LCOT  (I) 

Integer  which  indicates  whether  the  Ith  atom  has 
exited  the  sides  or  bottom  of  the  crystal,  has 
risen  from  the  crystal  surface,  or  is  still  within 
the  crystal. 


LD 

Number  of  the  next- to-hig'nest  numbered  mobile 
atom . 

LL 

Highest  numbered  atom  in  the  lattice. 

PA,PB 

Bohr  exponential  coefficients  in  Ar-Si 
interaction . 

PAC 

Ar-Si  exponential  force  coefficient  obtained  from 
dif ferentation  of  Born-Mayer  interaction 

potential. 


P AD, PBD 

Transition  potential  exponential  coefficients  in 
Ar-Si  interaction. 
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I 


PEXA/PEXB 

Born-flayer  exponential  coefficients  in  Ar-Si 
interaction . 


PFAD 

Ar-Si  Born-Mayer  force  exponential  coefficient 
obtained  from  differentiation  of  the  Born-flayer 
transition  potential. 

PFIV 

0.5 

PPPTC 

Ar-Si  force  evaluated  at  truncation  point  (ROEC)  . 

PFXA 

Ar-Si  force  exponential  coefficient  obtained  by 
dif ferentation  of  the  Bohr  potential  term. 

PPE(I) 

Potential  energy  of  the  Ith  atom. 

POT 

Potential  energy  generated  between  atoms  by  their 
interact  ion . 

PPTC 

Ar-Si  interaction  potential  evaluated  at 

nearest- ne ighbor  separation. 

PTC 

Si-Si  interaction  potential  evaluated  at 

near est- ne ighbor  separation. 
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mm 


PVBH 

Born-Mayer  term  of  the  Ar-Si  interaction 
potential. 

ROE 

Nearest- ne ighbor  distance. 

R0E2 

Nearest-neighbor  distance  squared. 

ROEA 

Cross-over  point  of  the  Bohr  and  Born-Mayer 
portions  of  the  Si-Si  interaction  approximation. 

ROEB 

Cross-over  point  of  the  Bohr  and  3orn-!1ayer 
portions  of  the  Ar-Si  interaction  potential 
approximation. 

ROEC 

Truncation  distance  of  the  Ar-Si  interaction 
potential  (set  at  nearest-neighbor  distance)  . 

R0EC2 

Ar-Si  truncation  distance  squared. 

ROEM 

Atomic  separation  distance  at  which  average, 
rather  than  instantaneous  force  between  atoms  is 
calculat  ed . 

SATE 

One-half  of  the  potential  energy  generated  by  the 
interaction  of  two  atoms. 
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RX(I)  ,RY(I)  ,RZ  (I) 

X,Y,Z  coordinates  of  the  Ith  atom. 


TPOT 

Total  potential  energy  of  all  atoms. 

XHAX,YMAX,ZMAX 

Maximum  X#Y,Z  coordinates  of  any  atom  which 
interact  with  the  atom  being  evaluated. 

XHIN , YMIN, ZMIN 

Minimum  X,Y,Z  coordinates  of  any  atom  which 
interact  with  the  atom  being  evaluated. 


could 


could 
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SIBP.OIFTNE  STEP 


THIS  SUBPOUTIN17  CALCULATES  THr  FORC”7 
T h I S IS  THE  nj  AM  2ND  L4TTICc  ADAPTION 
HIGH-SPEFO  VERSION  41TH  T HR E- -BRANCH 
OF  T H E FORCE  FUNCTIONS. 


9=TW=EN 

TF  th=  MUNDK'JR 
APDPPX:  MARIONS 


CCMM0N/COM1/°X(6OO) , RY(603) , R Z (600 ) , L CUT (600) , l_L,LD 
CCMM0N/C0M2/R0Et R0E2,RQFM,4C,PAC , PPTC , ptC , P CPTT , F DTC , 

IF* » PF  I V ,T  POT 

C0MMCN/C0M?/-X4 ,EXB , F XA , PEXA , Pr X B , PP XA , T Hc p * 

C OMMON / C 0 M5 / F X ( 600 ) , FY( 600) , FZ ( 6 00 ) , PP 7 ( 600 ) 

COMMON/ C0M22/RCF  A , R 0 3 , ROEC , R0  = 02 , CPO , C P L , C a2 , r P 3 , C c 0 , 
2CF1.CF2  »C  Gn  1 » C 30  2 , CGB 1 ,/‘GB2 
CCMMCN/C0M2A/PA, PB, PAO, PBO, A, Bf  PFACtCVED 

BOMBARD  I MG  ICN- LATTICE  ATOM  I NTcr  AC"  ! ""I  rOMP'J^'IDN 

1=1 

I F (LCUt ( I ) . NE . 3 ) GO  r0  203 

XMAX=RX ( ! ) +ROEC 

XMIN=RX ( I ) -ROEC 

Y M A X=  RY ( I ) +ROEC 

Y“ I N = RY ( I ) - ROEC 

ZMAX=RZ ( I ) +RCEC 

Z M I N = R Z ( T ) -ROEC 

DO  19  5 J=2,LL 

I F ( LCUT (J ) .Nc . 0 ) GO  '0  195 
IF (RY ( J ) .GE  .YMAX ) GO  195 
I F ( R Y ( J ).L=.YM!N)  GO  TO  195 
IF(RX(J J.GE.XMAX)  GO  T0  195 
I F(RX(J  ) .LE.XMI N)  GO  t3  ig5 
IF(PZ(J).GE.ZMAX)  Gn  T0  195 
I F ( R Z ( J ) .LE.ZMIN)  GO  tq  195 
QRX=RX( J ) — R X ( I ) 

DRY  =R Y ( J ) -R Y ( I ) 

DRZ=RZ( J)-RZ( T ) 

DI  ST2=DP  X**2+D5  Y**2  + ")P  Z**2 
IC(0IST2.GE  .R0EC2  ) GO  "0  195 
D I ST  = S0RT ( 0 IST2  ) 

IF(DIST.GT.ROEM)  GO  _0  165 
F 8M  = 6XP (PFXA  + PEX8*DIST  ) 

IF(OIST.GT.ROEB)  GO  t0  175 

fCRCE=CVcD*<1-03*DIST)*EXP(PA+PB*DISt)/DIST2+F3‘* 

GC  TO  180 

175  FORCE =E  XP ( PpA D+PBD*D I S T)  + C8M 
GC  tq  180 
165  DF  F=R  OFC-D I ST 

IF(DFF.lt.1 .3E-13  ) G°  T t 195 
FORCE  =(  EXP  ( PAC  + oEXP.*D  I ST  )-P  pptC  ) / DFF 
I F (FM.GT. FnRCE ) Gp  T0  195 
133  eCD  = FOR  CE/D!  ST 
F A =FOD- DP  X 
■ X ( I ) = F X ( I ) —FA 
F X ( J ) =F  X ( J ) +F A 
F A =FOD*DRY 
FY(  I ) =F  Y ( I ) -p A 
FY( J)=CY( J)+F4 
FA  =POD*DRZ 
FZ( I )=PZ( I ) -FA 
PZ(J)=FZ( J)+FA 
195  CONTTNUc 
233  CQNTINUE 

LATTICE  ATOM  INT-RACTION  COMPUTATION. 

DO  300  1=2, LD 

I F ( LGUT ( T ) .NE.3)  GO  ”0  333 

I P=T+1 

XMAX=RX ( I ) +ROE 
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XMIN=RX ( I )-ROE 

YMAX=RY( I )+ROE 

YM IM=RY ( I ) - ROE 

ZMAX=RZ(I > + ROE 

ZVIN=RZ ( I >-ROE 

DC  295  J=IP,LL 

I F ( LCUT  ( J ) • Nr  • 0 ) GO  T0  295 

IF (RY (J ) .GE  .YMAX ) GO  TO  295 

IF(RYU).LE.YMIN)  GO  T0  295 

I c ( RX ( J ) . GF  .XMA  X ) GO  TO  295 

IFIRXtJ  ).LF.XMIN>  GO  tq  295 

IF (PZ(J ) .GE.ZMAX)  GO  TO  295 

IF(RZ<J  l.LE.ZMIN)  GO  TO  295 

ORX=RX ( J ) -RX ( I) 

DRY=RY(  J ) - R Y ( I ) 

DRZ=RZ(  J)-P.Z(I) 

D !5.T2=DRX**2+DRY**2  + ORZ**2 
I F ( 01  ST 2 .GE  .ROE 2 ) GO  to  295 
DIST=S3R7 (DIST2 ) 

IF(OI ST .GT.ROEM)  GO  TO  265 
FORCE=  EXP ( FXA+EX3*0! ST ) 

IFIDIST.GT.ROEA ) GO  ”0  275 

FOR  Cc  =0 VE D~ ( 1-B*Q I ST) *EX® ( A + 3* 01 ST)/0!3T2  + C 
GO  TO  275 
2^5  DFF=RO=-DIST 

IF(DFF. LT.l.OF-10)  GO  To  295 
F0RCE  = (cXP(AC+FX5^DIST  )-PP TC ) /OFF 
IF(FM.GT.FORCc)  GO  to  295 
2 75  FCD=PORCE/OI$T 
FA=FOD*DRX 
FX( J)=FX( J)+PA 
FX(I)=FX( I ) -p A 
F A = FOO*QR  Y 
FY( J) =FY( J ) +p A 
F Y ( I ) =FY( I > -FA 
PA=FOD*DRZ 
FZ(J)=FZ( J)+FA 
FZ(I)=FZ(I) -FA 
295  CONTINUE 
3 J J CONTINUE 
RETURN 
ENC 
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Figure  13  - Ar-Si+4  Analytic  Potential  Function 
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(100)  Orientation , 1-keV  Sputtering  Events 
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Figure  25  - Sputtering  Mechanism  Trace 
[600  eV,  (110)  Orientation] 

121 

L 


ring  Mecha 
) Orientat 

7 


1 


j3®8®ffi§3®j 

^S3ffiL®~33 

®8®8®8®8 
&-&  - 'S'  -©— J 


8®8®8®8^ 

^8S8e8Sl 

^8SSS88 


g8R8_8Rg 

§3g3§=§=i 

bS8S8S8^ 

bS8S8S8® 


Figure  34  - (111)  Orientation,  600-eV  Sputtering  Events 
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Figure  36  - Sputtering  Mechanism  Trace 
[ 600  eV,  (111)  Orientation] 
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Figure  37  - (111)  Orientation,  1-keV  Sputtering  Events 
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Figure  42  - Copper  Single  Crystal  [ (100)  Orientation] 
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i Spread  in  Copper  (rimestep  26) 
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Timestep  3 


Figure  74  - 600 -eV  Precursor  Motion  in  Copper 
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Timestep  5 


Figure  75  - 600-eV  Precursor  Motion  in  Copper 
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Figure  76  - 600-eV  Precursor  Motion  in  Silicon 
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Timestep  5 


Figure  77  - 600-eV  Precursor  Motion  in  Silicon 
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